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Zur magnetischen Anisotropie von Einkristallen ae 
aus Eisen und Nickel 
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4 (Mit 16 Figuren) 


Das Problem der magnetischen Anisotropie von Einkristallen 
ist in den letzten Jahren mehrfach behandelt worden. Neben einer 
Zusammenfassung des bis zum Jahre 1932 vorliegenden Materials 
von R. Becker’) sind die Arbeiten von R. Gans zu erwähnen, die 
der theoretischen Deutung der von dem japanischen Forscherkreis ?) 
durchgeführten Versuche und der Beobachtungen von W.L. Webster?) 
und Sucksmith, Potter und Broadway*) gewidmet sind. 

Die genannten Überlegungen beziehen sich, ebenso wie die in 
der vorliegenden Arbeit durchgeführten, auf diejenigen Vorgänge, die 
bei hohen Feldern sich abspielen. Diese Vorgänge können dadurch 
beschrieben werden, daß man sich zu jedem Volumenelement einen 
Magnetisierungsvektor J, dessen Betrag gleich der Sättigungsmagne- 
tisierung J,, ist, zugeordnet denkt, der dann durch ein Feld von 
mehr oder weniger hohem Betrage gedreht werden kann. In dem 


Grenzfall unendlich großen Feldes liegt J in der Richtung des Feld- ‘= 
vektors H; mit abnehmendem Feld hat man im Rahmen dieser Bite 
Betrachtungsweise sich vorzustellen, daß der Vektor J sich all- ie 


mählich aus der Feldrichtung herausdreht in Richtung auf die 
nächstliegendste Kristallachse leichtester Magnetisierbarkeit, in der 
er ohne Feld liegen würde. Daß jedoch der genannte Mechanismus 
nicht mehr als rationelle Beschreibungsart verwendet werden kann R 
bei sehr kleinen Feldern, ist daraus zu ersehen, daß der Magneti- Nat 
sierungsvorgang dann in Form eines Barkhauseneffektes stattfindet, 
d. h. in Wandverschiebungen der homogen bis zur Sättigung magne- Br 
tisierten Bezirke’). Aus diesem Grunde ist ein Versagen der auf 
Drehungen des Magnetisierungsvektor beruhenden Beschreibungsweise 


1) R. Becker, Phys. Ztschr. 33. S. 905. 1932. 

2) K. Hondau.S. Kaya, Se. Rep. Tohoku Univ. 15.S. 721. 1926; K. Honda, 
H. Masumoto u. 8S. Kaya 17. 8. 111. 1928; S. Kaya 17. S. 639. 1928. 

3) W.L. Webster, Proc. Roy. Soc. (A) 107. S. 504. 1925. 


4) W. Sucksmith, H.H. Potter u. L. Broadway, Proc. Roy. Soe. (A) re ne 
117. 8. 476. 1928. 
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zu erwarten, wenn das Feld in der Größenordnung der Koerzitiv- 
kraft liegt. 

Von den oben erwähnten Arbeiten der Gansschen Schule ist 
in diesem Zusammenhang besonders hinzuweisen auf die beiden 
Arbeiten von R. Gans und Czerlinsky’); in ihnen erfolgt für Nickel, 
wo die Raumdiagonale Richtung leichtester Magnetisierbarkeit ist, 
die Behandlung der Magnetisierungskurve, d. h. die Berechnung der 
Parallelkomponente von J in Richtung von H für Würfelkante und 
Flächendiagonale, sowie die Behandlung der Normalkomponente von 
J in Würfel- und Dodekaederebene; beim Eisen, wo die Würfelkante 
die Richtung leichtester Magnetisierbarkeit ist, wird die Magnetisierungs- 
kurve in der Flächen- und Raumdiagonalen betrachtet, sowie die 
Normalkomponente von J in Würfel- und Dodekaederebene. Es sind 
also hiermit experimentell und theoretisch die einfachsten Fälle 
untersucht. In der Praxis erweist sich jedoch die Beherrschung 
dieser einfachen Fälle nicht als ausreichend. Andererseits genügt es, 
sich ebenfalls auf so hohe Felder zu beschränken, in denen der oben 
skizzierte Beschreibungsmechanismus zu Recht besteht. 

Die den folgenden Ausführungen zugrunde liegenden Versuche 
wurden an kreisrunden Scheiben ausgeführt, die an einem Torsions- 
draht mit ihrer Ebene parallel zur Feldrichtung aufgehängt waren; 
das von einem Elektromagneten erzeugte Magnetfeld selbst war um 
die Ebenennormale, d.h. die Richtung des Torsionsdrahtes, drehbar). 
Ist ©, der Vektor des Außenfeldes, so gilt allgemein, daß bei Nicht- 
parallelität von $, und J ein Drehmoment entsteht, das gegeben 


wird durch 
D= (dV 


wo [9,J] das vektorielle Produkt beider Vektoren bedeutet und 
über das Volumen der Scheibe zu integrieren ist. Bei der genannten 
Versuchsanordnung, bei der als Drehachse die Scheibennormale n 
festgelegt ist, kann nur die Komponente D des Drehmomentes in 
Richtung von n zur Auswirkung kommen; es ist dann bei homogener 
Magnetisierung und homogenem Außenfeld §,: 


| D=f{ J) aV = [J 
=f =H, [J,aV 


wobei J, die Normalkomponente von J in der Scheibenebene bedeutet. 


1) R. Gans u. E. Czerlinsky, Schrift. d. Königsberger Gel. Ges. Nat. Kl. 

9. 8. 1. 1932; Ann. d. Phys. [5] 16. S. 625. 1933. Zusammenfassung bei R. Gans, 
Leipziger Vorträge 1933, 91. 
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1. Allgemeine Grundlagen 

hh dem hier betrachteten Gebiet hoher Felder läßt sich die 
innere Energie eines Volumenelementes darstellen durch die Richtungs- 
kosinus «, der Winkel, die der Magnetisierungsvektor mit den Kristall- 
achsen bildet; sie läßt sich in erster Näherung darstellen durch die 
biquadratische Form 

Kia’? + + a”). 

Hierbei ist K eine Anisotropiekonstante, die bei Eisen positiv ist 
und etwa den Wert!) K =5-10° erg/cm® hat; bei Nickel ist sie 
negativ entsprechend der Tatsache, daß hier das Energieminimum 
in den Raumdiagonalen erreicht wird, und hat etwa den Wert K 
= — 5-10‘ erg/cm®. Physikalisch kann diese Energie gedeutet werden 
als Wechselwirkung zwischen dem Bahnmoment und dem Kreisel- 
moment des Elektrons). 

Es seien nun f, die Richtungskosinus des in dem Kristall vor- 
handenen Feldes 9; hierbei ist zu beachten, daß im allgemeinen 
diese Feldstärke von dem Außenfeld $, sich unterscheidet um ein 
entmagnetisierend wirkendes Zusatzfeld, das herrührt von dem in 
einem Aufpunkt P von der im PunktQ vorhandenen Magnetisierung Jo 
erzeugten Potential gp: 


wobei rpg den Abstand der Punkte P und Q bedeutet; das ent- 
magnetisierende Zusatzfeld hat dann den Betrag 


9, = — grad, pp = — [Jo gradg) grad, 
PQ 


d. h. 
1 
1,1) + fa V (J grad) grad — 


wobei jetzt jeweils nach den Koordinaten zu differenzieren ist, die 
als Integrationsvariable auftreten. Falls ein Zusammenhang von der 
Form J = x existiert, in dem x ein Skalar (Suszeptibilität) ist, so 
führt die letztere Beziehung durch Multiplikation mit x auf eine 
Integralgleichung für die Verteilung von g (W irschmidtsche Integral- 
gleichung). 

Herrscht in dem betrachteten Volumenelement die Feldstärke 9, 
so ist seine Energie 


129) USK (a,? + + — Joo (a, By B, +0; ß,)- 


Er R. Gans u. E. Czerlinsky, a.a. O. 
3 


2) F.Bloch u. G. Gentile, Ztsehr. f. Phys. 70. 8. 395. 1931. 
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Den Zusammenhang zwischen der Lage von 9 und der von J erhält 
man dann durch die Forderung, daß U als Funktion der «, zum 
Minimum zu machen ist unter der Nebenbedingung 

(1,3) = 

Setzt man zur Abkürzung 

(1,4) 


(1,5) — + @,?a@,? + @,? 
so lautet die Lösung 


as 


wobei :* der Winkel zwischen den Vektoren J und § ist; für ihn 
ergibt sich unter Benutzung der Beziehungen 


i=1 


die Darstellung 

(1,7) +yl-— 

wobei 

al A,?+A,?+ A,? = + a, 

— 4(a@,7 4,7 + 0,7 + a,” + a@,?a,*. 
Denkt man sich Polarkoordinaten eingeführt durch 

,=singcosp, a,=singsing, 4, = cose 

(1,8) | ß, =sino’cosgy’, 8, =sing’sing’, P,= cose 


wodurch die Nebenbedingung > Sa = 1 von selbst erfüllt ist, so 

stellt der durch formale Lösung des 
oU 

erhaltene Zusammenhang zwischen J und § dann das wirkliche 

Minimum dar, wenn 

(1,9) 4 


und entweder 


eu RU 
dg dgt 


oder 


iif 
We 
af 
= 
| 2 
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% 
= 
| 
= U= 7 (sin? 2 + sin‘ g sin’ 2 — HJ. 
(1,10) 
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und 
= K(2cos4o + sin?o (4cos?o — 1)sin?2 gm) + HJ. 

- (cos g cos + sin g sin 9’ cos — 9)), 
1,9) = 2Ksintocos4y + HJ. sin sin 0’ cos — p), 
= 2K sin’ cos sin 4g — HJ. cose sin sin(g’ — g). 


Die Achsen des Kristallsystems liegen nun im allgemeinen 
irgendwie zur Ebene der zu untersuchenden Einkristallscheiben orien- 
tiert; die Normale der Scheibe ist dann durch zwei Winkel o, 
und g, gegeben. Da bei der oben kurz beschriebenen Anordnung 
das Außenfeld $, in der Scheibenebene liegt, so kann die Lage 
von 9, beim Experimentieren leicht durch einen in der Blechebene 
gezählten Winkel o gekennzeichnet werden, der von der Projektion 
einer der Wiirfelkanten in die Blechebene ausgerechnet werden 
kann; er werde hierzu diejenige Würfelkante ausgewählt, von der 
aus der Winkel o gezählt wird; es bestimmen dann diese Würfel- 
kante, ihre Projektion in die Blechebene und 9, ein rechtwinklig 
sphärisches Dreieck, aus dem man hat: 


(1,8a) 9, = COS a SINO, 
(v9, und g, sind die Koordinaten von §,; 9, und g, die von n). 
Aus der Orthogonalität des Feldes und der Blechnormalen folgt: 
cos 9, COS + Sing, sing, cos(p, — y,) = 0; 
a. h 
(1,8b) cos(p, — P,) = — cotg o, cotg g,. 
Streng genommen handelt es sich um die Aufgabe, zu jedem 
Winkel o die zugehörige Verteilung der Magnetisierung in der 
Scheibe zu berechnen und daraus nach (1) die Größe des Dreh- 


momentes. Prinzipiell gesehen, könnte man sich hierzu den fol- 
genden Lösungsweg denken. Aus den Gleichgewichtsbedingungen 


ae =0, c= 0 folgen nach (1,6) und (1,7) die zwei Gleichungen: 
. C sin? 
coso’= cos#cosoe +CA,; sin(p’— yg) = — ing sin 4y 


oder symbolisch: 
= ®(o,y). 
Schreibt man eine Komponente der Beziehung (1,1) symbolisch in 
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so führt diese auf 
und mit (1,8a, b) auf 


wo das Symbol f o andeuten soll, daß o in dieser Gleichung auch a 


unter einem Integral vorkommt. Wegen (nO,)=0 hat man aus re 


(1,1) auch noch eine Beziehung: 


Die beiden Beziehungen G*=0, R= 0 a FR | 
gleichungen für die zu jedem Drehwinkel o gehörige Magneti- 
sierungsverteilung in der Scheibe. Man sieht hieraus, daß das 
Problem des Entmagnetisierungsfaktors und das der Theorie der 
Magnetisierungskurve ein Ganzes bildet und, streng genommen, nicht 
wie üblich, in zwei voneinander unabhängige Teilprobleme zerspalten __ 
werden kann. ” 

Um das hiermit skizzierte mathematische Problem zu umgehen, 
sei als erste Näherung angenommen, daß die Entmagnetisierung 
vernachlässigbar sei, also ®=,. Dann können als Zusammen- 
hang zwischen dem Drehwinkel o und den Winkeln o’ und gq’ ein- 
fach die Beziehungen (1,8a, b) angenommen werden, in denen 9, =’, 
Y,- y’ gesetzt werden. Die Magnetisierung der Scheibe wird dann 
als ortsunabhängig angenommen. 

Das Ergebnis des Versuches ist nun die Abhängigkeit des 
Drehmomentes D von dem Winkel o. Wir beweisen nunmehr den 
Satz, daß - - Stellen «, wo das Gleichgewicht gemäß (1,9) in- 
stabil wird, —— =o ist. Zum Beweis bilde man: 

dD de , öD d¢ 
ds da de 

Bezeichnet man nun mit dD/do’ eine Differentiation von D 
nach o’ nur an den Stellen seiner Darstellung, wo e’ explizit auf- 
tritt, so hat man: 

_ dD de , dD dg 
_ dD ög , 8D dg 
Aus den Gleichgewichtsbedingungen: 
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wobei natürlich gemeint ist: 


au 
folgt durch partielle Differentiation 
ö?U ög 6?U 


- +- 


Opdg dIpdg ag dg? dg dQ’ 


Hieraus erkennt man unmittelbar, daß dann, wenn 
wird, die Ableitungen 00/00’ und 0g/@o' unendlich groß werden; 
durch partielle Ableitungen der Gleichgewichtsbedingungen nach q’ 
beweist man dasselbe fiir 09/0’ und Oq/Oq’; also wachsen in 
dem Fall 4=0 auch 0D/00’ und OD/öyg' über alle Grenzen, 
womit zugleich die über dD/do behauptete Aussage bewiesen ist. — 
Wohlgemerkt ist bei dem letzteren Beweis keineswegs von der Dar- 
stellung (1,2) bzw. (1,10) der Kristallenergie Gebrauch gemacht; es 
ast nur vorausgesetzt, daß die Kristallenergie eine zweimal differenzier- 
bare Funktion ist. 


2. Das Drehmoment von Einkristallscheiben bei endlichem Feld 


Führt man Polarkoordinaten ein, gemäß (8), so erhält man für 
eine Scheibe, deren Würfelebene des Kristallgitters ¥" ar Scheibenebene 
6 
wegen g’ = 0 und damit = 0 sowie und 
die bekannte Parameterdarstellung: 


2HJ,, 


wobei die Magnetisierung der Scheibe als homogen angenommen 
ist. Die Elimination von o führt auf 


liegt, aus (1,10) und den Gleichgewichtsbedingungen © —=0, 


| sind o= 


- 2D 
40’ = arc sin TE + 4arc sin 


und unter Verwendung der Formel!) 


VHJ. 


oo 


4arc sin z= are sin (42 V1 — 2?) 1 — 422+ 424) 


1) Vgl. z.B. Madelung, Die math. Hilfsmittel des Physikers. 
1936, 3. Aufl., 8.54. 
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auf eine Darstellung von 49’ durch einen einzigen arc sin, in dessen 
D 


Argument das dimensionslos gemachte Drehmoment d= yg und 


die dimensionslos gemachte Feldstärke h = aoe eingehen; letztere 
Darstellung zeigt, daß das Drehmoment als Funktion von o’ die 
Periode 2/4 hat für jede Größe der Feldstärke. — Aus der zweiten 


Formel (2,1) folgt, daß man zu jedem o’ nur dann ein o erhält, wenn 


(2,2) h=0,5 


ist; diese Einschränkung liegt physikalisch daran, daß für sehr — 
kleine Felder keine Drehungen des Magnetisierungsvektors vor-- 
kommen. Eine allgemeine Gültigkeit der auf Drehung des Magneti- 

sierungsvektors gegründeten Theorie kommt demnach nur für Felder 


R 


360 


| 2 
0 90 180 270 —a’ 360 


Fig. 1. Gebiete des Minimums der Energie, wenn § und J 
in der Wiirfelebene liegen (K > 0) 


in Frage, die der Bedingung (2,2) geniigen; da andererseits zwischen 
dem Verhalten bei kleinen Feldern und dem bei groBen Feldern 
ein stetiger Übergang vorhanden sein muß, so ist zu beachten, daß 
eventuell unmittelbar oberhalb der durch (2,2) gegebenen Grenze 
der Drehprozeß noch nicht in Reinkultur realisiert zu sein braucht. 
Die Entscheidung, ob die betreffende Lage von $ und J ein wirk- 
| Minimum | der Energie liefert, wird wegen (1,9) gegeben durch: 
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mit weißem Röntgenlicht, in dem die Kobalt-K-Strahlung vorherrschte). 
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P=K|2cos4o + sin4o cotg (o’ — o)} >0 


sin 49 
sin (0° — 9) 


bo] bo] 


Q=K{2sinto + sing sin > 0. 


a Die Nullstellen und Pole von P und Q sind die von 


4tg(—o)=—tg4o und tg*o=tgo’. 


In Fig. 1 sind diese Kurven angegeben und die Gebiete, die 
bei Eisen keinem Minimum der Energie entsprechen, schraffiert. 

Fig. 2 zeigt die aus Formel (2,1) folgenden Kurven für das 
Drehmoment, wenn die Feldrichtung um die Scheibenachse gedreht 
wird. Man erkennt, daß im 
Gebiet 0,5 = h = 2 die Grenze 
des Gebietes, in dem U ein 
Minimum ist, erreicht wird 2p 94> 
und das Drehmoment der VA 06} 
Theorie nachspringen muß). 
Tatsächlich wurden solche 02 
Sprünge jedoch an einem Ein- 


kristall aus Hyperm 4, einer , 
Sonderqualitit des Krupp- 
schen Transformatorenbleches, 
nicht beobachtet, in Uberein- 46 
stimmung mit den Messungen 08+ 
von Honda und Kaya?) an -40Ll Ke 
_ Eisen. Zwar lag bei dem be- Fig. 2. 4 
nutzten Einkristall nicht die wenn (100) in Scheibensbene: a 
Wiirfelebene in der Blech- HJ, 
ebene, sondern die Ebene, 


deren Normale durch o, = 82,5°, 
y, = 10° gekennzeichnet ist (bestimmt durch Rückstrahlaufnahme 


Jedoch liefert die auf Drehungen des Magnetisierungsvektors 
gegriindete Theorie bei geniigend niedrigen Feldern auch dann 
Sprünge; so ergibt sich für die Dodekaederebene, wo g’ = 45° ist, 
unter Benutzung der Formel: 

sin 4g = 8sin g cos sin(45 + g)sin(45 — g) 

= 1) Vgl. auch die von N.S. Akulow in Wiss. Ber. d. Moskauer Staat 
universität H. 2. S. 8. 1934 gezeichneten Kurven. 


2) K. Honda u. S. Kaya, a. a. O. 
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durch ähnliche SchluBweise wie bei dem Fall der Würfelebene, 
daß auch g = 45° ist; dann erhält man: 


7 = sin — 0) 


und die bekannte Formel: 


Fig. 3. Gebiete des Minimums der Energie, wenn und J 
<n in der Dodekaederebene liegen (K > 0) 


Die das stabile Gebiet in Fig. 3 abgrenzenden Kurven ergeben sich 
gemäß (1,9) wieder aus P= 0, Q =0 und haben die Gleichungen: 
cotg(o’ — 0) = — te + sin 
bzw. 
sin 2 (1 3 sin’) 


cotg o’ = cotgo — 
tg ¢ 2 sin* o 


Fig. 4 zeigt die Veränderung des Drehmoments in Abhängigkeit von 
we _Feldrichtung; auch hier sieht man, daß in mittleren Feldern 
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Aus der Nichtexistenz der Sprünge ist zu schließen, daß eine 
allgemeine Gültigkeit der auf der Drehung des Magnetisierungs- 
vektors beruhenden Anschau- 


(2,5) h>2 10 


VK oat 
in Frage kommt. Für kleinere ae 


Felder ändert sich also dieDar- °%*7 
stellung (1,2) der Kristall- 04+ 
energie so, daß makroskopisch 2 
das Gleichgewicht während de 9 
Drehung des Feldes in dee 4, 


= 


Scheibenebene nicht labil wird }). 

, 
3. Die Magnetisierungskurve 0,67 

08+ 


Für Magnetisierung in 
einer Richtung, die der Wiirfel- = *9F 
ebene bzw. Dodekaederebene ~-%24 
angehört, ist die Gleichung Fig. 4. Drehmoment, wenn (110) 
der Magnetisierungskurve be- in Scheibenebene, Kurven h = const 
reits durch Formel (2,1) bzw. 

(2,6) gegeben?); die Kurven für den Fall des Eisens sind in den 
Figg. 5 und 6 dargestellt, in denen auch die aus den Sprung- 
bedingungen (2,3) und (2,4) sich ergebenden Koerzitivkräfte eingetragen 
| sind; aus der Bedingung (2,5) ist demnach nicht zu erwarten, daß 
die vorliegende Theorie im allgemeinen den Vorgang bereits für Felder 
in der Größenordnung der Koerzitivkraft und kleiner richtig darstellt. 

Bei Magnetisierung in beliebiger Richtung ist die Verwendung 
der oben angegebenen Formel (1,7) wenig geeignet, da sie zwar 
leicht aus gegebener Lage von J diejenige von § zu berechnen 
gestattet bei bekanntem Absolutbetrag von §; hier liegt jedoch der 
Fall gerade umgekehrt, daß die Lage von § bekannt ist und die 
von J gesucht. Zur Lösung der letzteren Aufgabe gehen wir aus 
von Formel (1,6); unter Benutzung der Werte der A, führt sie auf 


Ca a, B, — a, B, 


@, (a, — @,”) 


1) Uber einen Versuch, aus der hier beirachteten Beschreibungsart die 
Koerzitivkraft von Dauermagneten berechnen zu wollen, vgl. W.S. Messkin u. 
B. E. Somin, Zischr. f. Phys. 98. S. 610. 1936. 

2) Für den speziellen Fall der Flächen- und Raumdiagonale, vgl. bereits 
_N.8. Akulov, Ztschr. f. Phys. 67. S. 794. 1931; 69. S. 78. 1931; 81. S. 790. 1933; 
Gans u. E. Czerlinsky, a. a. O. 
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oder unter Benutzung der durch (1,8) angegebenen Polarkoordinaten: 


(3,1) C = sing’ ch Ba 


_ €08g sing’ — sin g cos@ cosg’ 


sin?g sing cospcos2q@ Sing cosg (cos?g—sin? g cos!g) ’ 


l 
1 0 1 


Magnetisierungskurven fiir verschiedene Richtungen 
der Wiirfelebene (K > 0) 


0 


N 
Fig. 6. Magnetisierungskurven für verschiedene Richtungen 
der Dodekaederebene (K > 0) 


woraus man die von H freie Gleichung 
, Cos 
COs @ Be 


1 — tg*o cos? 


, sin (g’ — 9) 3 
= 
4 tg 0 sin 4 n tg Q 


erhält, die sich mit y= cotgo umformen läßt in die bequemere Form 


sin + 39) sin 4 

2 


N 

| 

® 

+ 


H. Schlechtweg. Magnetische Anisotropie von Einkristallen usw. 585 


Bei bekannter Lage (0’, g’) von 9 läßt sich hieraus bei Annahme 
eines Wertes für g der Winkel o ermitteln und mit der aus (3,1) 
folgenden Formel 


der zugehörige Wert der Feldstärke. Da dies Verfahren sowohl bei EN 
Eisen als auch Nickel brauchbar ist, so ist bei der Lösung von (3,2) 
jeweils der stabile Zweig zu benutzen. Das Energieminimum liefert 
dabei derjenige Zweig der durch (3,2) definierten Funktion, der zu 
9 > liefert o—> 0’; hierbei geht dann h—> oo. Nach Gans’) 
ist h nicht aus der Feldstärke selbst zu berechnen, sondern aus der 
Differenz zwischen $ und derjenigen Feldstärke, die zur Erzeugung 
desselben J nötig ist bei Magnetisierung in der Richtung leichtester 
Magnetisierbarkeit. 


4. Das Drehmoment von Einkristallscheiben in sehr großen Feldern 


Da die allgemeinen Zusammenhänge, wie an dem Beispiel der 
Magnetisierungskurven in dem vorangegangenen Abschnitt gezeigt 
wurde, ziemlich verwickelt werden, betrachten wir die im Feld 
gedrehte Scheibe nur im Grenzfall sehr hohen Feldes, in denen eine 
Entwicklung nach negativen Potenzen des dimensionslos gemachten 
Feldes h angesetzt werden kann. Hier führt die Minimumforderung 
für den Ausdruck (1,10) auf?) 


Setzt man den die Kristallenergie darstellenden Anteil 


(sin? 20 + sintosin?2y)=f, 


ans u. E. Czerlinsky, a.a. 0. 


« 

/ 
| 

| 
| 

| 

1, A, & 

Einsetzen der Werte für die partiellen Ableitungen ergibt: 
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ae 


1 

: sin 4 — 5 sin 2 9’ sin 2 sin? 2 q’, 
= 
B= 1 4 ir 


A, = — A, (2 cos 49’ + sin?o’ (4 cos? 9’ — 1) sin? 2g’) 


sin*o’ sin 20’ 
+ 75 sin*e sin 20’ sin*4q’, 


B,= sin 20’ sin 4o’sin4g’ + sin? 20’ sin* 2 g’ sin 4q’ 


+ sin‘ 9’ sin 4g’ cos4q’. 
Zur Berechnung des Drehmomentes sind gemäß (1) die Komponenten 
des Vektors [OJ] zunächst festzustellen; sie sind 

= HJ... (sin 0’ sin p’ cos 9 — cos sing sin g), 
[HJ], = HJ. (cos sin g cos — sin 0’ cos g’ cos 
fant = [QJ], = H Je (sin sin sin(p — ¢’)). 
Mit Hilfe der Darstellungen 
cos 9 = cos(e — o’)cos — sin(g — sing’, 
sing = sin (0 — 9’) cos + cos(o — 0‘) sing’ 


der Entwicklungen 


1 


A, B, . 
y)= ‘sin — p) = 2+2, 


man 


(OI), =- 4, B, sin 0° cos 9’ cos +> 


x (94), =4, cos y’ — B, sin 0’ cos 9’ sin g’ ++ 


- sing’ — cos 9’ sin g’ (A, B,coso’+B, sing), 


Jh 


sin 0, cosp, + sin go, sin g, + 


| 
< 


ö = A, sing, sin(p,— 9) + B, cose, 
| if, Bo 
(4,3) 1 + h {( 4 sin 20 sın 0, sın ) 


sing’ 


er, 
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y 
(4,2) ı 
sin | A, 9 Sine cose Cosy (A, B, cos o’+ B, sing), 
age 
Dann wird 
aA 
Me 
oder : 
BE 
ER 
DE 


 kristallscheiben, deren 
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Diese Formel (4,3) gibt | a 
zusammen mit (4,2)und 4 
(1,8a, b) die Abhängig- 40% 
keit des Drehmomentes 
von dem Drehwinkel o 
des Einkristallbleches. 0 
Es ist ein beacht- 
liches Ergebnis, daß 
selbst für unendlich 
großes Magnetfeld das 
Drehmoment einen end- 
lichen, von Null ver- 
schiedenen Grenzwert 
annimmt, der durch 
das von 1/h freie Glied 
der Formel (4,3) gegeben 
wird. In den Figg. 7, 
8, 9 ist die Funktion 
lim ö(o) für einige kri- 
kh— oo 
stallographisch wichtige 
Ebenen, deren Miller- ay 
sche Indizes berechnet 
gemäß 


= +04 
jeweils angegeben sind, 4 
berechnet. Dem geo- I 
 metrischen Bild ent- 
nimmt man die Vermu- 
tung, daß sich an Em- % 


0080, = 


Ebene eine Würfelkante 
enthält, beim Heraus- 


&) drehen der F eldrichtung Figg. 7—9. Drehmoment 
Fu der Würfelkante hei verschiedener kristallographischer Orientierung 
dieselbe Neigung der der Scheibenebene 
Kurve (6) ergibt. 
Diese Vermutung läßt sich in der Tat beweisen; es folgt aus (4,3) 


ds dA, . „dp , dB 
= = 7, sing, sin (9,— — A, sing, cos(g,— + 6080, 


de ad 


\ \\ 
\ E 24 
/ \ \ 180° 
\ Wirftkane 
\ 
| 
Fig. 7 de 
4 
4 
BEN 
F \ EN) 
q 
\\ hf \ \ 
\ je \ \ j 
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und aus (4,2): 
‚d 
ae 


lo’ 

1 = 4cos4o’ + sın?2 sin? 2 


2 sin? cos 20’ sin? 2 + 2sin?o’ sin 2 0’ sin4g’ 


»b) folgt 


: Für die eine Ww ürfelkante 

o =o; fir a=0 wird 9 — om, = 180°; man erkennt a daB fiir 
solche Ebenen al. „endlich ist und ebenso an der Stelleo=0 auch 


— 21 = sin 2o’ sin 4 p + 4sin? 0’ cos 4 
do 
Es bleibt für die genannten Ebenen 
2 cos g’ cos 2 0’ cos Q, 
sin g’ 


dö„ COS Q, 


do 1 sin? o=0 


7 ein 08 (9 cog? sin?g, — 1) = 2(2 — 1) =2 
do o=0 YV1-— cos?eosin?o, 
Hiermit ist der Betrag der Neigung der 0, (o)-Kurve für die in Frage 
stehenden Ebenen ausgerechnet und seine Unabhängigkeit von der | 
Lage der übrigen Würfelkanten in Evidenz gesetzt. : 
Eine weitere, für die praktische Anwendung interessierende Größe 2 


nicht werden, in deren Ebene eine Wiirfelkante liegt: mit 9, = 90° 
o =o und — g, = 90°, d.h. 


erreic 


= 0 liefert fiir den Winkel o die Buben 


1 
-sin4o +37 1 sin 


_ 16-3sin?2 9. + V(16 — 3 sin* 2 g,)* — 8 (16 — 4 sin?2 Pn). 
2 (16 — 4 sin? 2 


ist in Fig. 10 dargestellt; man see, daß an den hier betrach- 


mi teten iam eine Bug: der Würfelebene um die in der Scheiben- 


“ad q 
iq 
a 
oder 
— 
Hi 
% 
— =. 
u 
£ 
Die Bedingung 
1" 
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ebene liegende Würfelkante den Wert des höchsten, der Würfelkante 
benachbarten Maximums nur wenig ändert (bis zu etwa 12°/,); das 


zweite Maximum ist das ya 


auf Orientierungsunter- es 
schiede schärfer an- 4 
2 >. ~ at 
An den Beispielen #2} 
Wy r a 
man stets eine viıer- 90 
zählige Symmetrie der Fig. 10. Höchstwerte des Drehmomentes En. a 


Kurven. Wann eine 
solche eintritt, erkennt 
man, wenn man die Potenzen des sin und cos in (4,2) ersetzt durch 
die sin und cos der ganzzahligen Vielfachen des Winkels. Zunächst 
sei bemerkt, daß, wenn o das Intervall von 0 bis a durchläuft, nach 
(Sa, b) die Variable q’ ein Intervall von der Länge a durchläuft, o’ 
jedoch ein kleineres Intervall. Man erhält 


in ihrer Ebene enthalten (0, = 90°, », beliebig) 


(4.4) 6, = ( sin 40’ + sin 20) sino, + gsin(4q' — 0’), 


sin 
wo g eine mit o variierende Amplitude und ö’ eine mit o variable 
Phase ist; 


2 sin? o sin? o (1 — cos 20’)? 
(= sin 40’ — — sin 20 cos? o. , 


1 sin4o’—2sin29’ sing 


tg 0 = 4 tg 0, . 


1 — cos 20’ sin 
Ist o, klein, so daß man sing, =~ 0 und cosg,= 1 setzen kann, 
d.h. praktisch eine Wiirfelebene in der Blechebene liegt, so sieht 
man aus (4,4) unmittelbar, daß man vier Nullpunkte von 0, (a) be- 
kommen muß; dasselbe erhält man fiir 9, = 90°; für alle übrigen 
Beispiele in Figg. 7, 8, 9 erhielten wir 4fache Symmetrie, was z. B. 
auch noch für die (321)-Ebene nachgeprüft wurde. Aus Stetigkeits- 
gründen ist daher anzunehmen, daß die Symmetrie von 0,, (a) stets 
4tach ist, mit Ausnahme der Oktaederebene, der gegenüber sich der 
Einkristall bei dem hier zugrunde gelegten Ansatz (1,2) praktisch 
so verhält, als ob er magnetisch isotrop wäre. 

Um die Verhältnisse an Eisen ausreichend zu beschreiben, nahm 
(Gans!) eine Erweiterung des Energieansatzes (1,2) vor, indem er 
die Kristallenergie in folgender Weise darstellte: 


(4,5) K (a@,* + @,? + Koa’w,’a,’. 
1) R. Gans u. E. Czerlinsky, a. a. O. fi) Ss 
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von Scheiben, die eine Würfelkante a 


— 
wee 2 
| 
j 
| 
a 
| 
| 
oz 
= Sr 
4 
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Die Wirkung dieses zusätzlichen Energieanteils ist am stärksten ie oe: 
der Gegend der Raumdiagonale; man erkennt dies besonders deut- 
lich an dem Schnitt der Energiefläche mit der Dodekaederebene, der in 
verschiedenen Maßstäben 
in den Figg. 11 und 12an- 
gegeben ist. In den durch 
(1,8) gegebenen Polarkoor- 
dinaten schreibt sich der 
zusätzliche Energieanteil 
ig Sin’ sin? 20 sin?2p. 
Beriicksichtigt man diese 
Korr mie nur in dem für 


0, 
180° 
Würfelkante 
Fig. 11. Schnitt der Fläche 4«,? a,* a,” 
mit der Dodekaederebene 


so ergibt sich, daß der mt 
ausgesprochene Satz über 
die Neigung der ö(o)-Kur- © 
0 ven in Ebenen, die eine 
| Würfelkante enthalten, 
| 

| 


o 


durch die Hinzufügung 
des Gliedes mit K’ nicht. 
geändert wird. 

Von der Wiedergabe 
der sich hiermit ergeben- 
den Korrektionskurven für 
das Drehmoment möge 
hier abgesehen werden. 


180° 

Würfelkante 5. Versuchsergebnisse 
Fig. 12. Schnitt der Fläche Der Vergleich zwi- 
+ 


Br schen Theorie und Ver- | 
mit der Dodekaederebene 


such wurde für einfachere 
Fälle der Würfel und Dodekaederebene bereits von Gans und 
 Czerlinsky ') durchgeführt. Da z. B. in der Wiirfelebene der 
Höchstwert des Drehmomentes unabhängig von H ist (Figg.2 und 4) 
gemäß der Theorie, so erstreckt sich der Gültigkeitsbereich der Theorie 
nur auf diejenigen Versuche, bei denen die von Hondo und Kaya”), 
ees und Sucksmith, Potter und Broadway‘) gemessene 

1) R. Gans u. E. Czerlinsky, a.a. O. 

2) K. Honda u. S. Kaya, a. a. O. 

3) W. L. Webster, a. a. O. 

4) W. Sucksmith, H. H. Potter u. L. Broadway, a.a. 0. 


. 
; 
— 
3 
\ 
If 
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Normalkomponente der Magnetisierung linear abhängt von .- In 


H, 
Fig. 13 sind die Versuchsergebnisse von Honda und Kaya an Eisen 
eingetragen, wobei die Stellen größten J, benutzt wurden. Man er- 
kennt, daß die Beobachtungen unmittelbar vor dem Sättigungsgebiet 
gemacht sind; denn der Steigungsfaktor der gestrichelten Geraden 
ergibt das zu unendlich großer Feldstärke gehörige Drehmoment, da 
sie in die aus den Meßpunkten gebildete Kurve praktisch ohne Knick 
einläuft. Die Versuche von Webster führen zu einer ähnlichen Dar- 
stellung. Aus Fig. 13 entnimmt man, daß demgemäß die Theorie 
etwa für das Gebiet H, > 1000 Oe gilt, d.h. rund h = H. “ve > 4, 

. 

/ “or 
j f Sucksmith 
/ / m) 
l N r 
00 7 0m 0 0,001 0,002 0.003 0000 0008 
Fig. 13. Auswertung der Ver- Fig. 14. Auswertung 
suchsergebnisse von Honda der Versuchsergebnisse an Nickel 


und Kaya an Eisen 


Für Nickel werden die Verhältnisse durch Fig. 14 wiedergegeben; 
die mit „Theorie“ bezeichnete Gerade ist gemäß Gans und Czer- 


linsky') berechnet mit X=5,12-10”* und — FE = 0,5. Die Versuche 


von Sucksmith, Potter und Ya ee sich für etwa 
H, > 2000, d.h. h, > 20 dieser Geraden; die von den Autoren an- 
gegebenen inneren Feldstärken wurden dabei mit Hilfe des Ent- 
magnetisierungsfaktors auf äußere Felder umgerechnet. Die Ergeb- 
nisse von Kaya?) liefern zwar ziemlich den richtigen Anstieg des 
geradlinigen Anfangsstückes, sind jedoch gegenüber der Theorie zu 


1) R. Gans u. E. Czerlinsky, 2.0. | 
29) K. Honda Kaya, aa. O. 
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kleineren Werten von J, verschoben; erhöht man J, um einen ent- 
H, >700 d.h. h, >7 bedeuten. Jedenfalls scheint sicher zu sein, 


gründete Theorie höchstens in dem ganzen Bereich von h, gilt, n 


nehmen zu können, und andererseits einiges Erfahrungsmaterial bei- 
zubringen über die Anisotropiekante K des mit 4,4 °/, Silegierten Trans- 
formatorenmaterials?). An Einkristallscheiben (20 mm Durchmesser, 
0,376 mm dick), die mit der Würfelebene des Kristallgitters in der Blech- 
ebene lagen, ergab, sich bei großem Außenfeld (1000—2400 Örsted) für 
die Abhängigkeit des Drehmomentes vom Drehwinkel o des Feldes die 
ar von der Theorie gemäß Fig. 2 geforderte Symmetrie. Da mit gegen oo 
; wachsendem Feld die Höchstwerte der Kurve D(c) einem Grenzwert 
zustreben müssen, so wurde dieser bestimmt durch Messung der 
Kurve Dyax/G in Abhängigkeit von 1/H, wobei G das Gewicht der 
Scheibe bedeutet. Fig. 15 zeigt das Versuchsergebnis (ausgezogene 
Kurve) und Fig. 15a die theoretisch, gemäß Formel (2,1) zu er- 
wartende Kurve (o’ = 22,5°); der Vergleich beider Figuren zeigt, daB. 
offenbar die versuchsmäßig aufgebrachten Felder groß genug waren, 
um die Richtung der in Fig. 15a gestrichelten Wendetangente in 
Fig. 15 noch experimentell zugänglich zu bekommen. Ist s das 

/ spezifische Gewicht (in unserem Fall gleich 7,6), so gilt für den 
u ee D von Dmax gemäß F ig. 15a und der Bedeutung von 0: 
= G 9818 

wobei s gemessen ist in g/cm? und K in erg/cm*; die Zahl D./G 
hat dann die Dimension cm. Man erhält somit 
K = 2,56 - 10° erg/cm?. 


Da dieser Wert erheblich kleiner ist als die in der Literatur bis 
Em: jetzt bekannten Werte von etwa 4,2 - 10° erg/cm?, ein kleiner Wert 
_ von K aber eine nur weiche Bindung der spontanen Magnetisierung 
an die Richtungen leichtester Magnetisierbarkeit bedeutet, so zeichnet 
sich das vorliegende Material durch eine besondere magnetische 


1) K. Honda u. S. Kaya, a. a. O. 
2) Sondermaterial „Hyperm 4“ der Fried. Krupp A.-G. mit besonders 


sprechenden Betrag, so würde dies eine Gültigkeit der Theorie fiir 
daß bei Nickel die auf Drehungen des Magnetisierungsvektors ge- % 


dem sie für Eisen gilt; die Grenze für das im Innern der Probe 
wirksame Feld wird wegen der Entmagnetisierung dann noch etwas __ 
herabgesetzt, so daß sie günstigenfalls wohl die Hälfte betragen mag, Br i 
wie man aus den Angaben von Honda und Kaya!) iiberschligt. 

Eigene Versuche dienten dazu, kritisch zu den Ansätzen dr 
Theorie und zur Art der Durchführung der Versuche Stellung 


ted d 
| 
>. 
tt wch 
es 
| 
7 
Pr 
: 
BR, großen Kristallen. 
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Weichheit aus. Für den Fall, daß nun Theorie und Versuch tat- 
sächlich von den gleichen Voraussetzungen ausgehen, darf der weitere 
Kurvenverlauf der Fig. 15 zu dem von Fig. 15a nicht in Wider- 
spruch stehen. Es wurde daher zu jedem Wert von Dy; der zu- 
gehörige Wert von 0 berechnet und aus Fig. 15 und 15a die be- 
treffenden Werte von 1/H und 1/k entnommen; da 4 = = + (Jo 
= 1550), so müßte sich aus H, h und K stets fiir J,, derselbe 
Wert ergeben; der untere Teil Fig. 15 zeigt, daB sich in keiner 


d 
053 


Too = 1550 Gauss 
K=26- 10° Ergfem® 


die Entmagnetisierung nicht vernachlassigbar ist im Gegensatz zu der 
im Verlauf von § 1 gemachten Annahme. Beriicksichtigt man sie 
angenähert, indem man von H den Betrag N J„ subtrahiert, wobei 
der Entmagnetisierungsfaktor N gemäß Snoek!) zu 0,267 sich be- 
rechnet, so erhält man für J in dem unteren Teil von Fig. 15 die 
aus H, = H — NJ,, berechnete Kurve (da in dem vorliegenden Fall 
für etwa Sr < 0,00146 stets h > 2 bleibt, so ist nach Fig. 5 hier 
das Gegenfeld NJ stets praktisch gleich N J.); man erkennt, daß 
bis etwa — = 2,10-* d.h. H > 5000, also etwa für h > 30 der 
berechnete Werte von J, sich als konstant ergeben würde. Nur 


1) J. L. Snoek, Physica 1. S. 649. 1934. 


= 

2 \ Einkristallschede 

ee) ; 511000 berectmet aus * \ 

\ 

Fig. 15 Fig. 15a 

a Weise Übereinstimmung ergibt, was offenbar davon herrührt, daß 7 2 

in 
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in diesem Gebiet ist also die übliche Art der rechnerischen Er- 
fassung der Entmagnetisierung statthaft, während in dem Bereich 


der experimentell nicht allzu schwierig zugänglichen Felder die n 
§ 1 skizzierte Rechnung in besserer Näherung durchgeführt werden __ 


müßte. Für die extrapolatorische Ermittlung von ö. scheint jedoch 
der entmagnetisierungsbedingte Unterschied zwischen der vorliegenden 
Versuchsdurchführung und der hier vorgenommenen Art der Aus- 
wertung der Theorie, der sich auch in der Lage des Wendepunktes _ 
und der Wendetangente bemerkbar machen muß, einstweilen noch 


nicht für die Praxis Fehler von einer Größe zu geben, die nicht — 


mehr tragbar wären. Legt man z. B. zwei Einkristallscheiben gleicher 
Abmessungen und er“ K, deren Ebene bei beiden 
die Würfelebene des Kristall- 

gitters ist, so erhält man 

torischen Grenzwert 0.., ob- 

wohl wegen des jetzt doppelt 
so großen Entmagnetisierungs- — 
faktors ) Wendepunkt und 
Wendetangente verschoben 

sind (vgl. die strichpunktierte 


Kurve von Fig. 15); stellt man — 2 


in diesem Fall aber 0 der als 
Funktion von 1/H,, so zeigen 
beide Kurven keine großen 
Abweichungen voneinander. 
An einem Einkristall 
gleicher Abmessungen wie 


die (123)-Ebene, d.h. eine 

Translationsebene?), lag, wurde 
die von der Theorie gefor- 

derte Symmetrie gefunden (vgl. 

Fig. 9); die für jeden der vier Extremwerte aufgenommene Abhängig- 

keit von der Feldstärke (Fig. 16) ergibt im Mittel K = 3,3 - 10° erg/cm?, 

also einen etwas größeren Wert als den Einkristallen der Würfel- 

ebene; bei der Berechnung von K wurde angenommen, daß der gerad- 


linig extrapolierte Wert gleich dem Grenzwert für unendlich hohes 
Feld sei, wodurch nach obigem K etwas zu groß ermittelt wird. 


1) J. L. Snoek, a.a. O. 
2) W. Fahrenhorst u. E. Schmid, Ztschr. f. Phys. 78. S. 383. 1932. 


oben, in dessen Blechebene 


- - 
> 
7 
r 
| 
i 


H. Schlechtweg. Magnetische Anisotropie von Einkristallen usw. 595 


Die aus K berechneten Grenzwerte für unendlich großes Feld (durch 
Kreuze in Fig. 16 bezeichnet), weichen höchstens um 9°/, ab gegen 
die extrapolierten Werte. Diese Abweichung hat ebenso wie wohl 
auch der höhere Wert von K seinen Grund mit darin, daß sowohl 
die geradlinige Extrapolation in Fig. 16 einen Fehler hereinbringt, 
als auch in der hier durchgeführten Auswertung nicht genau genug 
enthaltene Entmagnetisierung das Ergebnis fälscht; es wurde daher 
davon abgesehen, durch eine Korrektur nach Art von (4,5) die oben 
genannten 9°/, Abweichung zum Verschwinden zu bringen, da es 
zunächst noch fraglich erscheint, ob ihr für die hier durchgeführten 
Versuche eine physikalische Bedeutung zukommt. 

In welch hohem Grade ein Ansatz von der Form (4,5) ver- 
wendet werden kann zur Erklärung der auch vom Verfasser an 
zwei Einkristallscheiben beobachteten Tatsache, daß dann, wenn die 
Oktaederebene in der Blechebene liegt, sich die Scheibe im Gegen- 
satz zu der aus (4,3) für ö. zu ziehenden Folgerung im Magnetfeld 
nicht isotrop verhält, möge hier noch offen bleiben, jedenfalls könnte 
eine genauere Lösung des Entmagnetisierungsproblems zeigen, daß 
dann, wenn das Außenfeld in der Oktaederebene liegt, das im Inneren 
der Scheibe herrschende wahre Feld wegen der Kristallanisotropie 
nicht in der Oktaederebene liegt, womit zugleich, mindestens zu 
einem beachtlichen Teil, ein vermeintlicher Widerspruch zwischen 
dem Versuch und der auf den Ansatz (1,2) zusammen mit (1,1) ge- 
gründeten Theorie sich eventuell aufklären ließe. 


In einem ersten Abschnitt über grundsätzliche Fragen wird 
gezeigt, daß zwei Gruppen von Gleichungen zu lösen sind zur Aus- 
wertung magnetischer Versuche, die in dem Gebiet hoher Felder 
vorgenommen sind, wo der Vorgang durch eine Drehung des Magneti- 
sierungsvektors beschrieben werden kann; die eine Gruppe folgt aus 
dem Minimalproblem der Energie und liefert einen Zusammenhang 
zwischen dem inneren Feld und der Magnetisierung; die zweite 
Gruppe ergibt sich aus dem Problem der Entmagnetisierung und 
liefert einen Zusammenhang zwischen der räumlichen Verteilung von 
innerem Feld, Außenfeld und Magnetisierung. Eine exakte mathema- 
tische Behandlung steht noch aus. 

Es wird allgemein gezeigt, daß eine im Magnetfeld mit ihrer 
Ebene parallel zum Feld aufgehängten Einkristallscheibe ein sprung- 
haft sich änderndes Drehmoment erfährt an solchen Stellen, wo das 
Gleichgewicht instabil wird. An den einfachen Fällen, wo eine _ BIN: 
Würfel- oder Dodekaederebene in der Scheibenebene liegt, ergibt ager 


6. Zusammenfassung 


¥; 
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sich, daß Sprungstellen eintreten müßten, wenn der Magnetisierungs- _ 
vorgang bis zu sehr kleinen Feldern als Drehprozeß sich abspielen 
würde. Nach Angabe einer Formel für die Berechnung der Magneti- 
sierungskurve in beliebiger kristallographischer Richtung wird eine 
asymptotische Formel angegeben für das Drehmoment, das eine __ 
Scheibe in sehr großem Feld erleidet, und insbesondere der spezielle 
Fall von Scheiben betrachtet, in deren Ebene eine Würfelkante liegt. Es — 
werden Gründe angegeben, die dafür sprechen, daß der Mechanismus 

der Drehung des Magnetisierungsvektors im allgemeinen erst für 


Felder h = = Je > 4 gilt. Versuche zeigen, daß für das Krupp- 


sche Transformatorenmaterial die Anisotropiekonstante 
nur etwa 2,6 - 10° erg/cm® beträgt, also rg beachtlich weicher 
ist gegenüber dem z. B. von K. Honda u. S. Kaya benutzten 
schwedischen Eisen üblicher Verunreinigung. Es ergibt sich weiter- 
hin aus Unterschieden, die einerseits durch die Art der Durch- | 
führung der Versuche und andererseits durch die Art der Aus- | 
wertung der Theorie bedingt sind, von welcher Bedeutung das in 
der vorliegenden Arbeit nur in der iiblichen, fiir den vorliegenden 
Gegenstand etwas rohen, Weise beachtete Problem der Entmagneti- 
sierung ist. 
Essen, Physikal. Versuchsanstalt der Fried. Krupp A.-G., im 
1936. 
(Eingegangen 22. September 1936) te cat 
ene 


A 
ay 
2 
SER 
4 
3 
4 
a 


W.Denecke. Zur Methode der Sondenmessungen usw. 597 rs 


Zur Methode der Sondenmessungen 
in der Quecksilber- Niederdruck-Entladung*) 


Von Wollrath Denecke 


(Mit 16 Figuren) 


Inhaltsübersicht: A. Einleitung. Ziel der Untersuchungen. — 
B. Apparatur. — C. Meßergebnisse und ihre Diskussion: 1. Bestimmung der 
Geschwindigkeitsverteilung der Elektronen; 2. Bestimmung des Raumpotentials; 
3. Bestimmung der Elektronenstromdichten; 4. Bestimmung der Elektronen- 
konzentrationen; 5. Kontrolle der Wandladungsschichtdicken; 6. Verhalten 
der positiven Ionen. — D. Zusammenfassung. — E. Schrifttum. 


A. Ziel der 
Nach Einführung der dynamischen Sondenmethode von Lang- 
muir!=) ist eine große Anzahl von Arbeiten erschienen, in denen 
diese Methode angewendet und erweitert ist. Mit dieser Methode 
ließ sich der experimentelle Nachweis erbringen, daß die Elektronen 
in der positiven Säule einer Niederdruckentladung eine Geschwindig- 
keitsverteilung besitzen, die der Maxwellschen Verteilungsfunktion 
entspricht. Ein besonderes Ergebnis dieser neuen Untersuchungs- 
methode war es, daß unter gewissen noch nicht erfaßten Bedingungen 
im Plasma einer positiven Säule mehrere Geschwindigkeitsverteilungen 
der Elektronen nach Maxwell gefunden wurden. Dieses Ergebnis 
ist vielfach bestätigt worden. Seeliger und Hirchert®) fanden bei 
höheren Drucken, kleinen Stromstärken und langen Sonden in der 
Glimmentladung geknickte Charakteristiken, die sie durch Über- 
lagerung zweier Maxwellverteilungen der Elektronen zu erklären 
suchen. Mit kurzen Sonden erhalten sie stets nur eine Elektronen- 
temperatur, die jedoch höher lag als die mit langen Sonden ge- 
messene. Sie begründeten $-förmige Charakteristiken durch Wand- 
verunreinigungen, und können sie durch passende Behandlung der 
Gefäßwände vermeiden. Druyvesteyn’) hat durch 2 malige Differen- 
tiation der $-Charakteristiken, die im Niedervoltbogen in Ar und 
Ne gemessen sind, die Geschwindigkeitsverteilung der Elektronen 


*) Aus der von der Abteilung für Mathematik, Physik und Luftfahrt der 
Pe Hochschule Braunschweig angenommenen Dissertation. 
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bestimmt und kommt hierbei auf ein zweites Maximum der Maxwell- a 
schen Geschwindigkeitsverteilungskurve, dem er eine höhere Temperatur 2 
zuschreibt. Ähnliche Ergebnisse haben Sloane und Emeleus‘) 
erhalten, die auch vielfach geknickte Charakteristiken in Ar fanden. 

Bei langen Sonden erhalten sie zwei Temperaturen, bei kurzen da- — 
gegen nur eine und zwar die hohe. Sie führen dies auf die positive 
Ionenschicht zurück, die sich um isolierte Teile im Plasma befindet, 
wodurch erstens allgemein die gemessene Elektronenkonzentration im 
Plasma vermindert wird, und zweitens vorzugsweise die langsamsten 
Elektronen unterdrückt werden. Kömmnick®) findet im Hg-Nieder- 
druckbogen geknickte Charakteristiken, deren Reproduktion in be- 
stimmten Bereichen der Entladung jedoch zweifelhaft ist. Im all- 
gemeinen findet er nur eine hohe Temperatur. H. Paul?) findet 

in der geschichteten Säule in Wasserstoff im Schichtkopf ideale, n 
der Schichtmitte dagegen $-förmige Charakteristiken. Die Sonde 
wird hierbei durch die Anode in die Entladung eingeführt, sie wird 
durch eine lange Glaskapillare isoliert. Durch deren Aufladung. 
können aber Störungen in der Entladung eintreten. 

Infolge dieser Unsicherheit in der Bestimmung der eine’ 
und den daraus sich ergebenden Unsicherheiten in der Bestimmung — 
der Raumpotentiale, der Elektronen- und Ionenkonzentrationen er- 
schien es notwendig durch systematische Messungen auf experimen- 
tellem Wege die Fragen zu klären; vor allem die Frage: welche 
Geschwindigkeitsverteilung der Elektronen besteht wirklich in der 
freien ungestörten Entladung? Denn es erscheint unwahrscheinlich, 
daß bei Abwesenheit von störenden Wandladungseinflüssen in einem 
ee beschränktem Raume der 
positiven Säule gleichzeitig 

oder Elektronentemperaturen — 
bestehen können. 


B. Apparatur 
Die Messungen wurden 

in einer Hg- Niederdruck- 

entladung durchgeführt. Die 

Schaltung der Versuchsanord- 

Fig. 1. Schaltschema nung und die Ausbildung der 

und Anordnung der Sonde Sonde zeigt Fig. 1. Das seit- 

liche Hauptentladungsrohr, in 

ie die Sondenmessungen durchgeführt wurden, hatte einen Durch- 

messer von 5 cm und eine Länge von 80 em. Ds Katkodansmision 
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wurde bei allen Messungen konstant auf 5,0 Amp. gehalten. Durch Ver- 
änderung des Widerstandes zwischen Haupt- und Nebenanode (obere 
Elektrode) konnte die Stromdichte im Hauptentladungsrohr bei gleich- 
bleibendem Kathodenstrom geändert werden. Ein an der Kathoden- 
zuführung befestigtes W-Blech diente in bekannter Weise zur 
Fixierung des Kathodenfleckes. Die zylindrische Sonde aus W-Draht 
von verschiedenem Durchmesser ist an einem Stahldraht angeschweißt, 
der durch einen Barometerverschluß führt, so daß sie von außen her 
durch den Rohrquerschnitt verschoben werden konnte. Als Sonden- 
träger diente ein Glasröhrchen von einigen Zentimeter Länge und 
2,1 mm Durchmesser, das an seinem oberen Ende zu einer Kapillare 
ausgezogen war, durch die die Sonde gerade hindurchging. Am 
unteren Teil des Glasröhrchens war ein Eisenstück eingekittet; mit 
einem Magneten konnte somit der Sondenträger bewegt werden. 
Auf diese Weise gelang es, jede beliebige Sondenlänge einzustellen 
und gleichzeitig den Sondenträger an jede Stelle des Querschnitts 
zu bringen. In der bereits erwähnten Arbeit von Sloane und 
Emeleus°) wurde eine ähnliche verschiebbare Sonde benutzt, jedoch 
wurde diese mit einem Metallrohr umgeben, an das verschiedene, be- 
kannte Potentiale gelegt wurden. 

Während der Messungen lief dauernd die Diffusionspumpe. 
Der am Mc Leod-Manometer abgelesene Restgasdruck war kleiner als 
5.10% mm Hg. Die Kathode des Entladungsgefäßes tauchte in ein 
Wasserbad, das durch eine von Wasser durchflossene Kühlschlange 
dauernd auf konstanter Temperatur gehalten wurde. Durch Druck- 
luft wurde für gute Rührung und somit gleichmäßige Temperatur 
der Kühlflüssigkeit gesorgt. Der Hg-Dampfdruck betrug hierbei 
5.107? mm Hg. 

Der Sondenstrom und die Sondenspannung wurden mit Präzi- 
_ sionsinstrumenten gemessen, und zwar die Sondenspannung gegen 
- Anode. Bei gleichbleibender Stromstärke in der Hauptentladung 
wurde die Anodenspannung auf 0,1 Volt genau konstant gehalten, 
da es sich in den Vorversuchen herausgestellt hatte, daB schon 
kleine Schwankungen in der Anodenspannung zu großen Schwan- 


kungen in den MeBergebnissen führten. 


Es wurden drei grundsätzlich verschiedene Meßreihen quer zur 
Entladung bei jeweils gleichbleibendem Sondendurchmesser d auf- 
genommen. Bei allen drei Reihen waren die äußeren Betriebs- 


bedingungen konstant. 


1. Veränderung der Sondenlänge / und des Abstandes a des 
Sondenträgers und zwar so, daß die Sonde stets die obere Gefäß- 
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Der Sondenträger befindet sich in der Mitte des Entladungs- | 
rohres. Die Veränderung von / wird durch Herausschieben des. 
W-Drahtes bewirkt. D.h. a = const = 26 mm. 

3. Verschiebung des Sondenträgers bei unveränderter Sonden- 
länge | = const. 

; C. Meßergebnisse und ihre Diskussion 
1. Bestimmung der Geschwindigkeitsverteilung 
der Elektronen 

Die Bestimmung der Geschwindigkeitsverteilung der Elektronen 
in der positiven Säule der Niederdruckentladung geschieht aus der 
Neigung der im halblogarithmischen Maßstab aufgetragenen Charak- 
teristik Ini_=f(Vs). Der Elektronenstrom i_ wird aus dem ge- 
messenen Gesamtstrom i durch Subtraktion des positiven Ionen- 
stromes i, erhalten. Der bei einer negativ gegenüber dem Raum 
geladenen Sonde auf diese treffende Elektronstrom muß nach Lang- 
muir die Gleichung erfüllen: 


(1) i_=A-j_-e 
Hierin bedeutet: 
zur Sonde gelangender Elektronenstrom. 
- Sondenoberfläche. 
Elektronenstromdichte. 
Raumpotential. 
Sondenpotential. 
Elektronentemperatur, die der Maxwellschen Geschwindig- 
keitsverteilung der Elektronen entspricht. 

Bei der Ableitung dieser Formel ist von der Annahme aus- 
gegangen, daß die Geschwindigkeiten der Elektronen nach der 
Maxwellschen Verteilungsfunktion streuen. (Eine gute Übersicht 
über die wichtigsten Beziehungen, die bei den Sondenmessungen 
beobachtet werden müssen, findet sich z.B. im Handbuch der Physik!) 
und in dem Buch von Engel-Steenbeck"?). Die Gl. (1) gilt nur, 
so lange V; =V» ist. Unter dieser Voraussetzung besteht zwischen 
In:_ und Vz eine lineare Beziehung. Die Neigung der Geraden 
beträgt: 


8 
av,  K-T_ 
Hieraus die Elektronentemperatur T_ werden. 
Diese Elektronentemperatur kann durch die Beziehung 
7733 
in die mittlere Voltgeschwindigkeit umgerechnet werden. 
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In der Literatur sind häufig Beobachtungen angegeben, wonach 
die halblogarithmischen MeBergebnisse Ini_ = f(Vs) nicht mehr einer 
Geraden entsprechen*—'*), Es zeigen sich Knicke in dem Verlauf, 
woraus man zunächst rein formal auf 2 Temperaturen bzw. Elek- 
tronengeschwindigkeiten geschlossen hat, wobei aber auf die Ent- 
stehung zweier Maxwellverteilungen im Einzelnen nicht weiter ein- 
gegangen wurde. Da aber anzunehmen ist, daß in einem beschränkten 
Raum der positiven Säule in genügender Entfernung von den Ge- 
fäßwänden nur eine Elektronengeschwindigkeit bestehen kann, soll 
in Folgendem versucht werden, festzustellen, welche der bei den 
Messungen auftretenden Elektronengeschwindigkeiten wirklich im 
Plasma vorhanden ist. 

Bei vielen der hier angestellten Versuche traten ebenfalls Ab- 
weichungen von der idealen Form der Charakteristik auf. Der 
Strom ö_ stieg bei einem bestimmten gegenüber dem Raumpotential 
noch negativen Potential viel stärker an, als es nach dem bis- 
herigen Verlauf zu erwarten war. Es ergaben sich also zwei ver- 
schiedene Neigungen der Charakteristik, aus denen nach der obigen 
Beziehung zwei verschiedene Temperaturen ermittelt werden konnten. 
Diese Abweichungen von der idealen Form der Charakteristik waren 
bei gleichbleibenden Versuchsbedingungen sehr gut reproduzierbar. 

Bei Sondenlängen bis ungefähr 3 mm waren keine Knickpunkte 
in der &eraden festzustellen. Diese traten erst bei größeren Längen 
hervor. Der Sondendurchmesser war hierbei 0,8 mm, der Durch- 
messer des Sondenträgers 2,1 mm. Bei einem anderen Sonden- 
durchmesser von d = 0,2 mm traten die Knickpunkte bei ungefähr 
denselben Längen auf. Der Durchmesser des Sondenträgers war 

in allen Fällen derselbe. Nach den Arbeiten von Sloane und 
Emeleus°) war zu schließen, daß die auf dem Sondenträger befind- 
lichen Ladungen an dem Auftreten der zwei Temperaturen beteiligt 
sind. Diesen Einfluß der Wandladungen zeigt folgender Versuch. 
; Der Sondentriger befindet sich in der Mitte der Entladung 
- (a= 26mm). Die Sonde selbst (d = 0,8 mm) wird 5 mm lang ge- 
macht, so daß gerade deutlich die beiden Geschwindigkeiten auftreten. 


Hierauf wird die ganze Anordnung, also Sonde und Sondenträger, 
soweit auf die Wand zu verschoben, daß die Sonde die obere 


Gefäßwand berührt (a = 6,5 mm). Es zeigt sich nur wieder eine, 
und zwar die große Geschwindigkeit. Fig. 2 zeigt die MeBpunkte. 
- Um das Auftreten der 2 Temperaturen auch in der letzten Stellung 
zu erreichen, muß die Sondenlänge ungefähr 7 mm betragen. Dieser 
Versuch zeigt, daß in der Nähe von Wandladungen aus den Messungen 
2 Geschwindigkeiten abgeleitet werden können. 
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Die Gefäßwände laden sich, wie jeder in die Entladungsbahn — 
gebrachte isolierte Körper, gegenüber der eigentlichen Entladung _ 
negativ auf. Es bilden sich daher hier Raumladungsschichten aus, | 
deren Dicke später berechnet 
pi werden soll. Hier werden für die 
} er weiteren Betrachtungen nur die aus 


d=0,8mm / den Messungen sich ergebenden 
Werte zugrunde gelegt. Aus dem 
Kurvenverlauf der Raumpotentiale 
in Fig. 13 (vgl. später) ist zu er- 
pee f sehen, daß bei einer Sondenlänge 
ER / von ungefähr 7 mm bei d = 0,8 mm 
Aa H der Einfluß der Wand nachläßt. 
RG if (Bei d = 0,2 mm wird dieser Ein- 
4 i, fluB erst bei 1 = 8 gering). Der 
Stérbereich der Wandladungen 
beträgt also 7+8 mm. Diese 
Schicht besteht sowohl an der 
/ Gefäßwand wie auch an dem 
of / | Sondenträger. 
Fig. 2. Sondencharakteristiken halb dieser so bestimmten Schicht 
bei gleichem / und verschiedenem a eine Störung der Geschwindig- 
keitsverteilung der Elektronen ein- 
treten kann, kann man nach dem Beispiel von A.v. Engel und 
M. Steenbeck (Bd. 2, 8. 23) die Relaxationsstrecke S berechnen. 


Nach Engel-Steenbeck folgt: 
Swl 108. - at 


Fiir zwei verschiedene ie ergeben sich die Werte 
der Tab. 1. 


Relaxationsstrecke S und Schichtdicke x 
T N x S 
mm °K iw cm cm 


0,2 | 9700 2,2 0,8 0,8 
0,8 1,2 0,7 1,3 


Für die Messungen bei d = 0,2 mm ist die errechnete Relaxations- 
strecke gleich der Schichtdicke x, während sie bei d = 0,8 mm 
gréBer ist. Nach den Angaben von Engel-Steenbeck kann sich 
aber eine Maxwellsche a unter den 
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tronen ausbilden kann. 


sitzen. 


erreichen können. 
innerhalb der Schicht läßt sich nicht 
angeben. Als obere Grenze muB 
die Geschwindigkeit gelten, die 
durch die Größe der Wandaufladung 
(im vorliegenden Falle 5 Volt) ge- 
geben ist. 
Auf eine solche aufgeladene 
Wand können also nur diejenigen 
Elektronen gelangen, deren kine- 
ist, wenn U der Potentialabfall 
innerhalb der Schicht ist. Wird in 
diese Schicht eine Sonde gebracht, 
so können offenbar nicht mehr alle 
Elektronen diese erreichen; die 
He langsamsten werden durch das Feld 


Je näher die Sonde 
gebracht wird, um so 
Elektronen können sie 


rückgetrieben. 
der Wand 

weniger 

erreichen. 

verteilung 


- 101! em?, 
wurde. 


in der Schicht abgebremst und zu- 
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Elektronen schon auf einer Weglänge ausbilden, die nur 50°/, des 
errechneten Wertes beträgt. Unter dieser Annahme würde in beiden 
Fällen die Schichtdicke x größer sein als S/2, so daß sich inner- 
halb der Schicht eine neue Maxwellverteilung unter den Elek- 
In die Schicht können alle 
eindringen, die eine dazu normale Geschwindigkeitskomponente be- 
Infolge der Feldwirkung von der Wand her, wird aber den 
Elektronen ein Teil ihrer kinetischen Energie entzogen, so daß die 
langsameren vollständig abgebremst wurden. 
Geschwindigkeit innerhalb dieser Schicht wird also sicher größer 
sein als der in der freien Entladung, da ja nur die schnelleren 
durch die Schicht hindurchkommen und die Sonde bzw. die Wand 
Ein genauer Wert für diese Geschwindigkeit 


Elektronen 


Der Mittelwert der 


7.=9000 °K 
N.=21*10" cm? 
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verteilung der Elektronen 


so wird diese in 


Wenn also auch in der Entladung eine Maxwell- 
der Elektronen besteht, 
zs solcher Wandladungen nicht mehr voll bestehen. 


der Nähe 
Die Fig. 3 


wie sie aus gemessenen Werten bei d= 0,2 mm 
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Die mittlere Geschwindigkeit der Elektronen liegt bei 1,2 Volt. 
Die Wandaufladung beträgt, wie wir später noch sehen, ungefähr 
5 Volt [J.v.Issendorf!?) nimmt sogar 10—15 Volt an]. Die Wand 
selber können nur die Elektronen erreichen, deren Geschwindigkeit 
größer als 5 Volt ist. Wie man aus Fig. 3 sieht, ist dieses nur ein 
sehr kleiner Teil sämtlicher Elektronen. 
Aus diesen Betrachtungen folgt, daß sich bei kleinen Sonden- 
längen die gemessene Geschwindigkeit der Elektronen mit Annäherung 
an einen isolierten Körper stark ändert; denn je kleiner die Sonde 
wird, bzw. je mehr sie sich der Gefäßwand nähert, um so stärker 
wird der Einfluß der Wandladung auf die Sonde, d.h. um so mehr 
wird der Schnitt in Fig. 3 nach rechts nach größerer Geschwindig- 
keit verschoben. 
Auch durch eine Potentialflächendarstellung läßt sich die 
scheinbare Zunahme der Elektronengeschwindigkeit mit abnehmender 
Sondenlänge erklären (vgl. Fig. 4). 
ze Die Potentialflächen sind für den 
triger Fall gezeichnet, daß die Sonde sich — 
gerade auf Raumpotential befindet. 
Die Sondenlänge ist kleiner als die 
sich um den Sondenträger ausbildende 
Schichtdicke. Es muß jetzt eine Po- 
tentialfläche geben, die sich in der 
Nähe der Sonde teilt. Diese Teilung _ 
tritt bei um so negativerem Potential — 
auf, je kürzer die Sonde ist und | 
0-1 -2-3-4-5 -§ -4#-3 -2-1 0 umgekehrt. Die Potentialfliche, die 
* 4. Schematische Darstellung diese Teilung zeigt, gibt auch F ; 
der Potentiallinien. V, = V,  Yade das Potential an, gegen das 


8° die Elektronen anlaufen müssen, 
auf die Sonde zu gelangen. Man erkennt aus dieser Darstellung, ey: 
da die Potentialschwelle um so größer ist, je kleiner die Sonde 


wird. Es können also bei sehr kleinen Sonden nur die schnellsten 
Elektronen die Sonde erreichen, die langsamen werden durch die 
davor gelagerte Potentialschwelle von der Sonde ferngehalten. Die 
ülektrunenstromdichte wird also bei verschiedenen Sondenlängen 
immer verschieden sein (solange die ganze Sonde noch innerhalb 
a Schicht ist), und zwar wird sich die Stromdichte bei ver- 
schiedenen Sondenlängen entsprechend der in der freien Entladung 
herrschenden Geschwindigkeitsverteilung der Elektronen ändern, — 
da bei sich ändernder Sondenlinge die Potentialschwelle ver- 
schieden ist. 
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Mit diesen Überlegungen stimmten die Versuchsergebnisse 
durchaus überein (vgl. Fig. 8. Bei Annäherung der Sonde an die 
Gefäßwand (bei gleicher Länge / = 1 mm) wächst das Raumpotential 
stark an, außerdem wurde die Neigung der halblogarithmischen 
(Geraden immer geringer; d.h. die Elektronengeschwindigkeit wird 
größer. Gleichzeitig sinkt bei gleicher Sondenspannung der auf die 
Sonde treffende Elektronenstrom; d.h. verringert sich die Elektronen- 
konzentration. 

Nehmen wir jetzt eine Sonde, deren Länge / größer ist als die 
Dicke x der Wandladungsschicht. Wir können uns dann diese 
Sonde in 2 Teile zerlegt denken. Der eine Teil befindet sich inner- 
halb der Schicht, während der andere sich in der von der Wand- 
ladung des Sondenträgers unbeeinflußten freien Entladung befindet. 
Zur Vereinfachung machen wir auch noch die Annahme, daß inner- 
halb der Schicht überall die gleiche Elektronenkonzentration und 
Geschwindigkeitsverteilung herrscht. 

In der freien Entladung haben die Elektronen eine Geschwindig- 
keitsverteilung entsprechend einer Temperatur 7T,, in der Schicht 
herrsche die Temperatur T,. Entsprechend werden die Längen, 
Stromdichten und Ströme bezeichnet. Der 
gesamte auf die Sonde fließende Strom 
besteht aus der Summe der beiden Einzel- |, 
ströme. Für Fall 2 (vgl. S. 600, a = const | (--—- 


= 26mm) ergibt sich das Bild der Fig. 5. 


eV 


| 
~ 


I 
a =j,-l,-n-d-e = | 
ig, 5 


/ eV eV 

i=a-d (1, “ji .e ER +L-j, 
Für den Gesamtstrom i ergibt sich aus dieser Betrachtungsweise 
eine Summe von 2e-Funktionen. Ist 1 = /,, d.h. gleich der Schicht- 
dicke x, so tritt nur die Temperatur T, in Erscheinung. Mit wach- 
sendem ! kommt immer mehr der Einfluß von T, mit hinzu. Es 
ist zu erwarten, daß sich der Knickpunkt mit größer werdendem 1 
nach links verschiebt, da jetzt der Einfluß von !, größer wird. Die 
Tabellen zeigen die Ergebnisse für d=8mm. Die erste Spalte 
gibt die einzelnen Sondenlängen an. In der zweiten sind die Raum- 
potentiale angegeben, die sich mit wachsender Söndenlänge ändern 
müssen. Die dritte Reihe gibt an, bei welcher Sondenspannung V, 
der Knickpunkt auftrat. In der vierten Reihe endlich sind die 
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a Tabelle 2 Tabelle 3 
Abhängigkeit des Knickpunktes V, 
der Differenz V,—V, von der Sondenlänge 
‘Fall Fall 2 
a und | variabel a = const, | variabel 


Ve Ve |%-% 
Volt | Volt Volt 


4,90 

4,20 
9 4,10 
13 4,05 
185 3,95 
245 | 3,90 


wo 


7 
‚so 


or or 
or or 


16 13 
17% 

1 


= Im Fall 1 (vgl. S. 600) haben wir nicht x, sondern 2x einzusetzen, — 
My er hierbei ständig zwei Wandladungen wirksam sind. Wie zu 
erwarten, tritt bei dieser Anordnung der Knickpunkt erst bei größeren 
Sondenlängen auf im Fall 2 der größten Länge im Falle 2 


b) 
Fig. 6. Auftreten und Verschiebung des Knickpunktes bei verschiedenem ! | 


 (l= 24,5 mm) muß er wieder nach rechts wandern, weil jetzt die Sonde 
Jin die Wandladungsschicht der Gefäßwand hineinragt. Bei diesen 
Überlegungen ist zu berücksichtigen, daß die Sonde nur ein mittleres 
Z 5 umpotential anzeigt. Längs der Sonde ist ja sicher ein konstantes 
Potential, während in der Aufladung längs dieser Strecke ein solches 
nicht immer bestehen wird. Sloane und Emeleus®) haben bi 
ihren Untersuchungen mit einer Doppelsonde ähnliche Ergebnisse 
erhalten. Als Sonde benutzten sie einen Molybdändraht, der durch 
einen Cu-Schutzring hin und her bewegt werden konnte. Andiesen _ 
Schutzring konnten verschiedene Spannungen gelegt werden. Wurde _ 
er bei kleiner Sondenlänge (1 = 1mm) mit der Anode verbunden, © 
so trat in der Charakteristik der Knickpunkt auf. Wurde unter — 
sonst an ihn ein stark 
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gelegt, so konnte aus der Charakteristik nur die große Geschwindig- 
keit bestimmt werden. Die Erklärung hierfür ist folgende: Solange 
der Schutzring mit der Anode verbunden ist, können dort keine Ionen- 
raumladungen entstehen, die Elektronen können ohne Geschwindig- 
keitsverlust ungehindert die Mo-Sonde erreichen. Bei stark negativem 
Potential an dem Schutzring entsteht dort aber eine solche Schicht, 
durch die die langsamsten Elektronen von der Sonde ferngehalten 
werden. Der Schutzring entspricht in diesem Falle — bei negativem 
Potential — dem Sondenträger der vorliegenden Arbeit. Die Er- 
gebnisse dieser Arbeit bestätigen also vollkommen die ähnlichen 
Untersuchungen von Sloane und Emeleus. Darüber hinaus soll 
aber auch ein Verfahren angegeben werden, wie auch bei kleinen oder 
verschiedenen Sondenlängen die wirklich in der Entladung vorhandene 
Elektronentemperatur bestimmt werden kann. 

Eine einfache Überlegung, die gleichzeitig die Grundlage der 
Langmuirschen Sondenmethode ist, führt hier zum Ziel. 

Wir bezeichnen mit o_ die mittlere ungerichtete Geschwindigkeit 
der Elektronen im Plasma; j_ sei die Dichte des hierdurch hervor- 
gerufenen Elektronenstroms; ist außerdem n_ die Konzentration der 
Elektronen, so gilt: 
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(2) j_=14-..n_-V_ 


Die mittlere ungerichtete Geschwindigkeit ö der Elektronen bleibt 

im Plasma an allen Raumstellen dieselbe. Nach dem Boltzmann- 

schen Satze stehen die Elektronenkonzentrationen n, und n, an 

2Raumstellen 1 und 2, deren potentielle Energien sich um den Betrag 
eV 


&- V unterscheiden, im Verhiltnise *-?. Da nach Gl.(2) die Strom- 
dichten j_ den Konzentrationen n_ proportional sind, so wird die 
Stromdichte j an einer Raumstelle 1, deren potentielle Energie um 
e-V gegenüber einer anderen Raumstelle 0 mit der Stromdichte 7, 
kleiner ist 


Legen wir für die Messungen eine ideale Sonde zugrunde, d. h. eine 
punktförmige Sonde ohne Isolation, dann müssen im halblogarith- 
mischen Maßstabe die Stromdichten für die dazugehörigen Raum- 
potentiale auf einer Geraden liegen, die der Gl. (3) entspricht. Aus 
der Neigung dieser Geraden kann eine Temperatur T_ berechnet 
werden, die der wahren Geschwindigkeitsverteilung der Elektronen 
in der Entladung entsprechen muß, denn die einzelnen Raumpotentiale 
sind ja nichts weiter als Stellen mit verschiedener potentieller 
Energie. Die einzelnen Charakteristiken müssen in diesem idealen 
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= 


Falle zusammenfallen, man erhält so den in Fig. 7 dargestellten 
Verlauf. Die Charakteristiken unterscheiden sich nur durch das 
Raumpotential, bei dem das Sättigungsgebiet beginnt, und somit auch 
SA durch den Wert der Stromdichte. 
Messen wir mit einer wirk- 
lichen Sonde an mehreren Stellen 
des Querschnitts, so erhalten wir 
4 für die einzelnen Raumpoten- 
tiale V, verschiedene Stromdich- 
ten j_. Bei konstanter Sonden- 
länge / sind die Stromdichten den 
Strömen direkt proportional. Es 
zeigt sich nun, daß durch die ein- 
zelnen Raumpotentialströme sehr 
gut wieder eine Gerade gelegt 
werden kann, deren Neigung die niedere Temperatur angibt. 

: (Vgl. hierzu Figg. 8 und 9.) Diese Gerade entspricht durch- 


= 


ce» 


~ 


wachsender 
Abstand a 


ds 
Fig. 7. Charakteristiken 
einer idealen Sonde 
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Fig. 8. Fig. 9. 

Bestimmung der wirklichen Geschwindigkeitsverteilung der Elektronen 

ae im Plasma bei konstanter Sondenlinge. i_ = /(V,) 


aus der Gl. (8), aus ihrer Neigung kann also die in der freien 
Entladung herrschende Geschwindigkeitsverteilung der Elektronen 
bestimmt werden. Aus Fig. 8 folgt für d = 0,8 mm, T = 9330° K 
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= 1,217. Aus Fig.9 folgt für d= 0,4mm, T = 9300° K = 1,20V. 
Besonders erwähnenswert ist hierbei die Tatsache, daß in den einzelnen 
_ Charakteristiken diese niedere Temperatur gar nicht in die Erschei- 
nung tritt, da die Sonde sehr klein gehalten wurde. (J = 1 mm). 
Die Tabellen 4 und 5 geben die Zahlenwerte für die Figg. 8 
and 9. Zu beachten ist der starke Anstieg des Raumpotentials in 
der Wandnähe und der damit zusammenhängende Abfall der Elektronen- 
stromdichten. 

3 Tabelle 4 Tabelle 5 
d=08mm, / = 1 mm d = 0,4 mm, / = 1 mm 
Elektronentemperatur, Raumpotential und Stromdichte als Funktion 
des Abstandes der Sonde von der Gefäßwand bei gleicher Länge / 


= Vi j- xe) j- 

°K Volt | Amp./em~? mm | °K Volt | Amp./em— 
25 31200 | 5,5 0,022 sey 2,5 | 37700 | 6,8 0,005 
5.3600 4,4 0,083 7 = 4 24900 5,5 0,025 
| 21900 | 4,0 0140 6,5) 24200 | 5,15 0,039 
26 | 21100 | 3,6 0213 11 | 22900 | 47 0,074 
A 26 22600 | 4,4 0,109 


stand a des Sondenträgers ergeben “” 7 F 
sich jeweils andere Raumpotentiale J-, /: 
(vgl. weiter unten) Man bestimmt 7” see fl 

jetzt fiir mehrere Sondenlingen das A 
jeweilige Raumpotential und die dazu A 

gehörige Stromdichte. Diese Strom- 
dichten werden im halblogarithmi- 

schen Maßstab über den entsprechen- “” +7 
den Raumpotentialen aufgetragen. / 

4 A Die einzelnen Werte müssen wieder 

der Gl.(3) genügen. Es ergibt sich, 

= wie in Fig. 10 gezeigt ist, wieder eine 2 / 
fear aus deren Neigung abermals / 

die kleine Geschwindigkeit bestimmt — 

wind. Fig. 10 Verhältnisse “8 ~6 


Fig. 10. Bestimmung der wirk- 
lichen Geschwindigkeitsvertei- 
lung der Elektronen; im 

der 
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tronen übereinstimmend angegeben wird. Tab. 6 gibt für ver- 
schiedene Sondendurchmesser die nach diesem Verfahren bestimmten 
Elektronentemperaturen, wobei bei den einzelnen Meßreihen die 


Sondenlänge verändert wurde. 


Tabelle 6 


Elektronentemperaturen bei verschiedenen Sondendurchmessern 


Arten bestimmen. Bei langen Sonden erhält man oberhalb des 
Knickpunktes schon den richtigen Wert, bei kurzen Sonden muß 
man das oben angegebene Verfahren benutzen. Ein Vergleich zwischen 


den derart bestimmten Elektronengeschwindigkeiten ergibt, daß die ® 
einzelnen Abweichungen untereinander auch bei mehreren Sonden- 5 


/ 


-0 -5 


Fig. 11. Bestimmung der 
wirklichen Geschwindig- 
keitsverteilung der Elektro- 
nen im Plasma bei verschie- 
denen Hauptstromstärken. 


durchmessern nicht größer sind als 10°/,. 

Bei veränderter Bogenstromdichte 
muß sich an derselben Stelle der Ent- 
ladung ein anderes Raumpotential ein- _ 
stellen, d.h. an derselben Stelle haben 
wir jeweils eine andere potentielle Energie. 
Bestimmen wir also mit einer Sonde an 
derselben Stelle der Entladung bei ver- 
schiedenen Hauptstromdichten die Raum- 
potentiale und die hierbei auftretenden 
Stromdichten, so müssen diese beiden 
Werte in demselben Zusammenhang 
stehen, wie oben für verschiedene Ab- 
stände angegeben war. Es muß also 
auch hier gelten: 

Fig. 11 zeigt die Verhältnisse für die 
Sondenabmessungen : 
d=0,8 mm, /=8 mm. 

Die Verbindungslinie der Raumpoten- 
tialstréme ergibt wieder die niedere 
Temperatur. (T_ = 9200° K = 1,19 Volt.) 
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Aus dieser Figur ist noch zu erkennen, daß der Einflußbereich der 
Wandladungen bei verschiedenen Bogenstromdichten verschieden ist. 
Bei kleinen Hauptstromdichten ist er größer als bei großen. Daher 
tritt bei kleinen Stromdichten nur die große Geschwindigkeit auf, 
während bei großen Stromdichten der Knick in der Charakteristik 
sichtbar wird. Dieselben Messungen sind auch bei kleinen Sonden- 
längen durchgeführt. Das Ergebnis ist genau das gleiche, nur daß 
hierbei in den einzelnen Charakteristiken nur die hohe Geschwindig- 
keit auftritt. Die Neigung der Verbindungslinie der Raumpotential- 
ströme liefert auch hier wieder die kleine Elektronengeschwindigkeit, 
die mit der nach dem Verfahren der Figg. 8—10 bestimmten Ge- 
schwindigkeit übereinstimmt. 

Auf S. 603 hatten wir die Schichtdicke x an der Gefäßwand zu 
ungefähr 7 mm angegeben. Die Stromstärke im Hauptladungsrohr 
betrug hierbei 3,0 Amp. Wird diese auf 1,0 Amp. herabgesetzt, so 
muß die Schichtdicke entsprechend ansteigen. In einer Annäherungs- 
rechnung kann man sagen: 

j, = positive Ionenstromdichte. Sie ist direkt proportional der Haupt- 
bogenstromdichte J. 7, ~ J. Also können wir setzen: 


I~}, m 


Bei 1,0 Amp. Bogenstromstärke wird also die Schichtdicke um den 
Wert Y3 größer als bei 3,0 Amp. zı4= 12 mm. In diesem Falle 
(!= 8 mm) ist die Länge ! kleiner als die Schichtdicke x, es kann 2. 
also nur eine und zwar die große Elektronengeschwindigkeit auftreten, BT 


wie auch die Messungen richtig ergaben. Re 
2. Bestimmung des Raumpotentials 
Das Raumpotential Vin der positiven Säule einer Niederdruck- 


entladung bestimmt man am einfachsten aus der im halblogarith- Pe: 
mischen Maßstab aufgetragenen Charakteristik Ini_=f(Vs). Slangeg 

Vs < Vz ist, ergibt sich eine Gerade. Sobald aber V;,=V, wird, >i 
treten hiervon Abweichungen auf, wir erhalten ein Sättigungsgebiet. 
Der zur Sonde gelangende Strom i_ muß nun nach den verschiedenen 
Raumladungsgleichungen berechnet werden. Das Potential, bei dem 
diese Abweichung auftritt, ist das gesuchte Raumpotential. Zu be- 
merken ist, daß die Sonde stets nur einen mittleren Wert des längs 
ihrer gesamten Länge herrschenden Raumpotentials anzeigt. Es 
genügt daher nicht allein die Angabe der Sondenlänge, sondern es 
muß auch der genaue Stand der Sonde angegeben werden. Bei dem 
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Auftreten von mehreren Geschwindigkeiten ist zur Bestimmung des 
Raumpotentials stets der oberste Knickpunkt zu wählen. Diese Art 
der Bestimmung des Raumpotentials ist unabhängig von der Gestalt 
der Sonde; sie kann also in jedem Falle angewandt werden. Für 
Kugel- und Zylindersonden gibt es noch eine zweite Bestimmungs- 
möglichkeit, die aus den Raumladungsgleichungen folgt. 

FR Für Zylindersonden gilt: 


4= 2 .— . 
(1 + ) 


* oe - 


> 


K-T_ 
Extrapoliert man aus dieser linearen Beziehung mischen ı und V, 
auf den Werti— = 0, so erhält man: i 


2 


“sea 


Bei Kenntnis der Elektronentemperatur T_ kann man auf diese Weise 


das Raumpotential bestimmen. 
Für Kugelsonden folgt ähnlich: 


Auch hier erhält man wieder bei Kenntnis T_ für i_ = 0 das Raum- 
potential. 

Bei den Untersuchungen wurde eine Zylindersonde benutzt, 
deren Länge / verlängert werden konnte; für 1>d haben wir es mit 
einer wirklichen Zylindersonde zu tun; wird / = d, so kann die Sonde 
angenähert als Kugel betrachtet werden; wenn 1<d, so liegt sehr 
angenähert der Fall der ebenen Sonde vor. Wie weit die einzelnen 
Bereiche gelten, wird an Hand der Messungen dadurch bestimmt, 
daß zunächst rein formal für jede Meßreihe die Bestimmung des 
Raumpotentials nach den Formeln für die ebene, Kugel- und Zylinder- 
sonde durchgeführt wurde. Es wurde also so verfahren, als ob die 
Sonde wirklich eine Ebene, ein Zylinder oder eine Kugel war. Hier- 
bei bestätigten sich durchaus die oben angestellten Überlegungen. 

Bei den Rechnungen für die Kugel- bzw. Zylindersonde geht 
in die Gleichung, aus der das Raumpotential bestimmt wird, die 
Elektronentemperatur mit ein. Hier muß natürlich die in der Ent- 
ladung wirklich vorhandene, also niedere Temperatur eingesetzt 
werden. Fig. 12 zeigt die Verhältnisse bei einem Sondendurch- 
messer d = 0,8 mm. Das Raumpotential Vz ist in Abhängigkeit von 
der Sondenlänge | aufgetragen. Die Anordnung war so, daß die 
Sonde stets die obere Gefäßwand berührte und nur der Sondenträger 
nach unten verschoben wurde, d.h. die Länge der Sonde stetig zu- 
nahm (Fall Die Auswertung der Meßreiben unabhängig 
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von der Sondenlänge nach allen 3 Methoden durchgeführt. Zur 
Übersicht wurden auch noch die Werte K- T_/e eingetragen und 
zwar sowohl für die hohe wie für die niedere Temperatur. Bei 
großen Längen | zeigen die Bestimmungen des Knickpunktes der 
halblogarithmischen Geraden gute Übereinstimmung mit den Er- 
gebnissen der Berech- 5 
nungen nach der Zylin- 4 
dersondenmethode. Für 4! 
kleinere Längen bleiben 
diese Rechnungswerte 
hinter den Meßergebnis- 
sen zurück, statt dessen 
haben wir Ubereinstim- ; 
mung der Messungen 
und der Rechnungswerte 

nach der Kugelsonden- 
methode. Für sehrkleine 9g 5 70 6 
Längen haben auch diese 
Rechnungen zu kleine 
Werte und wir können uns nur noch auf die Bestimmung des 
Kniekpunktes der In i_ — V,-Geraden verlassen. Diese Methode 
liefert bei jeder Länge die richtigen Ergebnisse. Das wahre Raum- 
potential ist durch die stark ausgezogene Kurve gegeben, die mit 
der für die ebene Sonde erhaltenen Kurve zusammenfällt. Es ist 
aus den gemessenen Werten möglich, den Knickpunkt in der halb- 
logarithmischen Geraden mit einer größten Abweichung von 0,1 Volt 
zu finden. Bei der Berechnnng nach der Zylindersondenmethode 
können Fehler namentlich durch die Temperaturbestimmung ein- 
treten. Die höchsten Abweichungen hierbei betrugen aber 10°/,, 
so daß bei einer Elektronentemperatur von T_ = 95000° K ein 
Fehler von weniger als 0,1 Volt für die Bestimmung des Raum- 
potentials besteht. Die Abweichungen zwischen diesen beiden Me- 
thoden bei langen Sonden liegen innerhalb dieser Fehlergrenzen. 
Aus Fig. 12 folgt also, daß bei d= 0,8 mm von 8 mm an die 
Berechnung des Raumpotentials nach der Zylindersondenmethode zu 
richtigen Ergebnissen führt. Für kleinere Längen ist nur die Be- 
stimmung aus dem Knick der In 1_— V,-Geraden zuverlässig. 

Für die in der Beschreibung der Apparatur angegebenen drei 
Fälle wurden die Raumpotentialmessungen bei einem Sondendurch- 
messer d= 0,8 mm durchgeführt. Im Fall 1 muß sich bei kleinen 
Sondenlängen ein größeres Raumpotential ergeben als im Fall 2, 
weil bei ersterem 2 negative Potentiale auf die Sonde einwirken, 
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” Kuge/- ” 
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während im 2. Falle nur die Wandladung des Sondenträgers Ein- 
fluB haben kann. Fig. 13 zeigt die erhaltenen Ergebnisse. Im | 
Fall 3 betrug die Sondenlänge 1 mm; als Abszisse ist jetzt nicht , 
sondern der Abstand a zu nehmen. Aus den Kurven J und 2 kann 
der Einfluß der Gefäßwand erhalten werden, wenn man berück- 
sichtigt, daß in beiden 
% | d=08mm Fällen der Einfluß des 
| Sondenträgers auf die 
m === Sonde für gleiches / der- 
selbe ist. Die Differenz- 
2 bildung 1—2 schaltet die 
Wirkung des Sonden- 
trigers aus und gibt | 
0 5 00 —2, “a) 2% jetzt nur den Wirkungs- 
Fig. 13. Einfluß der Gefäßwand. bereich der Gefäßwand = 
Ve = fi) bzw. f(a) an. Mit dieser Differenz- 
bildung darf man jedoch 
nicht so weit gehen, daß bei Fall 2 die Sonde bei großem / schon 
in den Bereich der Gefäßwandladungen kommt. Bei ungefähr 
7 mm ist der Einfluß der Gefäßwand verschwunden; daß hierbei die 
Differenzkurve noch nicht den Wert 0 erreicht hat, ist dadurch zu 
erklären, daß der Sondenträger im Fall 1 weiter in die Entladung 
hineinragt, wodurch sicher die ganze Entladung gestört wurde. 
Diese Störungen machten sich bei Hauptstromstärken von 0,6 Amp. 
(0,033 Amp./cm”?) so bemerkbar, daß die Anodenspannung etwas 
anstieg, wenn der Sondenträger bis an die obere Gefäßwand ge- 
bracht wurde. 

Die Kurve für den Fall 3 (= const= 1 mm) muß im all- 
gemeinen für alle Werte außer /= 1mm höher steigen als die 
beiden anderen Kurven, da ja bei jeder Stellung die negative Auf- 
ladung des Sondenträgers auf die ganze Sonde einwirkte. Die Werte 
für a wurden hierbei insofern korrigiert, als a den Abstand von 
der inneren Gefäßwand bedeutet, um so zu dem gleichen Anfangs- 
punkt für die Kurven 7 und 3 bei / = 1 mm zu kommen. 

Bei den ganzen Messungen ist zu bemerken, daß die Sonde 
stets nur einen Mittelwert des längs ihrer ganzen Länge herrschenden 
Raumpotentials anzeigt; in diesen Mittelwert geht natürlich auch 
immer der Einfluß des Sondenträgers mit ein. Je länger die Sonde 
gemacht wird, um so weniger kommt dieser Einfluß zur Geltung. 
Es ist günstig, den Durchmesser des Sondenträgers sehr klein, da- 
gegen den der Sonde möglichst groß zu wählen. Dieser Sonden- 
durchmesser ist aber durch die Bedingung begrenzt, daß der Ent- 
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ladung nicht zu viele Träger entzogen werden dürfen, da sonst 

hierdurch wieder eine Störungsquelle geschaffen wird. ate er 

3. Bestimmung der Elektronenstromdichten 

Die Bestimmung der Stromdichten geschieht im Punkte des 
Raumpotentials. Aus Gl. (1) folgt, wenn Vs = Vp ist, 

= = A 

Es ist also der zur Sonde gelangende Elektronenstrom i_ durch die 

Sondenoberfliche A zu dividieren, um j_ zu erhalten. Wenn wir 

mit einer idealen Sonde in einer von keiner Wandladung beeinflußten 

Entladung messen, so wird die Stromdichte j_ für alle Längen / der 


Sonde dieselbe sein; da sich aber um den Sondenträger, sowohl wie 
um die Gefäßwand Wandladun- 


gen ausbilden, so ändert sich, 
wie wir sahen, mit wachsender 
Länge das Raumpotential (mit 
wachsendem / wird es kleiner). N 
Damit wird aber nach der in 
dem Kapitel Temperaturbestim- 1 & 
mung näher abgeleiteten Be- et 
ziehung a ar + 
+ ca 
j_ größer. Bleibt das Raum- 4 
potential bei verschiedenen 
Längen konstant, so bleibt auch 7? 0 
j_ konstant. Im Fall 1 fiel das 


Raumpotential ständig ab, folg- Fig 14. 
lich muß ein ständiger Anstieg 
von j_ stattfinden. Fig. 14 zeigt die Abhängigkeit der Elektronenstrom- 
dichte von der Sondenlänge für d = 0,8 mm. Die aus den einzelnen 
Messungen folgenden Werte sind durch Kreuze gekennzeichnet. Die 
Abweichungen sind auf ungenaue Bestimmung des Raumpotentials 
und vor allem der Sondenfläche zurückzuführen. Die eingezeichnete 
Kurve ist aus der Fig. 10 und 13 konstruiert. 


4. Bestimmung der Elektronenkonzentration 

Bei Kenntnis der Geschwindigkeitsverteilung der Elektronen 
kann die Konzentration berechnet werden. 
_ /2a-m 13 
n_ = = 4,03 - 10 VT 
ist hierbei in einzusetzen. 


ae vag 
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Weiter oben wurde darauf hingewiesen, daß beim scheinbaren 
Auftreten von zwei Geschwindigkeiten nur die kleinere in der Ent- 
ladung wirklich besteht. Diese ist daher bei der Konzentrations- 
bestimmung einzusetzen. Der Verlauf der Konzentrationen in Ab- 
hängigkeit von der Sondenlänge (Fig. 15) ist genau der gleiche wie 
der der Stromdichten, so daß 
hierauf nicht noch einmal ein- 
gegangen wird. Es ist bei 

+ kleinen Sondenlängen nur dar- 
BEE; auf zu achten, daß auch wirk- 
lich der richtige Temperaturwert 
eingesetzt wird. Da beikleinen _ 
Sondenlängen in der Charakte- 

wee ristik nur die große Geschwin- 

re digkeit auftritt, kann hier leicht 
die Konzentration zu niedrig 
ermittelt werden; es ist daher 
bei den Auswertungen stets die 
mm niedere Geschwindigkeit ein- 

Fig. 15. n_ = fi gesetzt worden. Die Kreuze 
in Fig. 15 geben wieder die 

einzelnen Meßergebnisse an, während die ausgezogenen Kurven 
durch Rechnung aus Fig. 15 nach obiger Beziehung erhalten 
wurden. Zu beachten ist der starke Abfall der Konzentration in 
der Nähe der Wand. Da mit der Sondenmethode in unmittelbarer 
Nähe der Wand und an der Wand selber keine brauchbaren Er- 


die Verhältnisse an der Gefäßwand durch eine Rechnung erfaßt — 
werden. 


5. Kontrolle der Wandladungsschichtdicken ies JE 


Die Verhältnisse an der Gefäßwand können durch die Sonden- 
methode, wie sie zu den vorliegenden Messungen benutzt wurde, 
nicht erfaßt werden. Wir sind hierbei auf Rechnungen angewiesen, 
die nach den Angaben von Engel-Steenbeck'*t durchgeführt 
werden sollen; es sollen hierbei unter Benutzung der Meßergebnisse 
zahlenmäßige Angaben über die Konzentrationen, Potentiale usw. 
an der Wand gemacht werden, um so zum Schluß die Schicht- 
dicke x dieser Wandstörungen errechnen zu können. Die vor den 
Gleichungen in Klammern () angegebenen Zahlen beziehen sich auf 
die fortlaufende Nummerierung im Engel-Steenbeck, Bd. 2, Die 
a im n Engel- Deuabnck sind unter der Voraussetzung 


gebnisse mehr erzielt werden können, sollen im folgenden Kapitel _ 
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abgeleitet, daß die freie Weglänge klein ist gegenüber den Rohr- 
dimensionen. Im vorliegenden Fall ist aber diese Voraussetzung Bar 
nicht erfüllt; es muß also an Hand der MeBergebnisse zunächst —_ 
geprüft werden, ob noch Übereinstimmung der Rechnung mit den 
Meßergebnissen vorliegt. Die Figg. 13 und 15 zeigen die Abhängig- 
keit der Raumpotentiale bzw. der Konzentrationen von der Sonden- 
länge. Hier besteht die Beziehung: 

K-T_ 

€ N, 

Die Differenz V,—V, ergibt sich aus Fig. 13. Bei kon- 
stantem V, (für ] = 24,5 mm) betrug N, = 7,2. 10!° mm”?,. Für 
verschiedene Werte von V, wurde jetzt N, berechnet. Die Er- 
gebnisse sind in Fig. 15 als Punkte eingetragen. Man sieht, daß 
eine hinreichende Genauigkeit zwischen Rechnung und Messung 
vorliegt. Es ist also berechtigt, die Rechnungen auf den vor- 
liegenden Fall anzuwenden. 

Aus den Gl. (160), (157), (154), (170), (179) und (171) von 
Engel-Steenbeck (Bd. 2) läßt sich die Elektronenkonzentration 
an der Wand bestimmen zu: 


Nr= 2,45 - 10° cm”?, 
- Nach Gl. (181f) folgt somit für das Wandpotential: 
= Va— V, = — 493 Volt. 


0 


(181¢) 


- 


ar 


Da die Annahme gemacht war, daß jeder Träger, der die Wand Ts 
erreicht, dem Gasraum endgiiltig entzogen ist, daB also die Anlage- en 
rungswahrscheinlichkeit = 1 vorliegt, so kann aus der Wandkonzen- === 
tration der zur Wand gerichtete Strom berechnet werden. Diese = 


Stromdichte beträgt: 
= 6,4. 10-* Amp./em-?.; 


> 


Nunmehr kann nach den bekannten Raumladungsgleichungen 
für zylindrische Anordnung aus je und Vr— V, die Dicke der 
Raumladungsschicht an der Gefäßwand berechnet 


wenn jr in Amp./cm?, V in Volt und r in cm eingesetzt FE 
r= 1,7 ( cm. 

Der Rohrradius betrug 2,4 cm; der Radius der Raumladungs- ee 


Somit wird die Schichtdicke selbst: x = 0,7 cm. Ea x 


= 
> 
% 
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Dieser Wert stimmt sehr gut mit dem schon vorher aus den Mes- 
sungen sich ergebenden Wert überein. 


6. Verhalten der positiven Ionen 


Bei stark negativem Sondenpotential wird der Sondenstrom 
nur durch positive Ionen getragen. Wir erhalten in diesem Bereich 
eine lineare Beziehung zwischen Strom und Spannung. Für kleinere 
negative Spannungen überlagert sich diesem Ionenstrom der Elek- — 
tronenstrom. Der Ionenstrom im Raumpotential kann in Annähe- © 
rung durch Extrapolation erhalten 
werden, wie Fig. 16 zeigt. Für 
die „Temperatur“ der positiven 

Ionen lassen sich mit Sicherheit — 
nur zwei Grenzen angeben. Die — 
obere Grenze ist die Temperatur _ 
ets der Elektronen, die die Ionen je- — 
SD doch infolge ihrer zu geringen Be- | 
weglichkeit nicht erreichen werden; _ 
die untere Grenze ist die Tem- 
peratur des neutralen Gases. In 
Annäherung nimmt man für die 
re. T, der Ionen die halbe Elektronentemperatur an. 


Fig. 16. Bestimmung des positiven 
Ionenstromes bei V,, 


kleiner Sonde wird der Unterschied noch stärker, so daß also n i 
der positiven Säule eine starke positive Raumladung vorhanden sein 
müßte. In anderen Arbeiten®®) treten sogar noch größere Unter- 
schiede auf. Da derartige Raumladungen in einem hochionisierten 
Gase nicht möglich sind, so bestünde einmal die Möglichkeit, daß der 
positive Ionenstrom durch die Sondenmethode nicht richtig gemessen 
wird. Hierzu ist zu sagen, daß der genaue Abfall des positiven 
Ionenstroms in der Nähe des Raumpotentials durch den mehrere 
Zehnerpotenzen größeren Elektronenstrom nicht bekannt ist; es 
wird aber unwahrscheinlich sein, daß durch diese Fehlbestimmung 
ein derartiger Unterschied in der Konzentration hervorgebracht 
wird. Die in der Literatur angegebenen stärkeren Unterschiede 
zwischen n_ und n, können zum Teil dadurch begründet werden, 
daß zur Bestimmung der Elektronenkonzentration die hohe Tem- 
peratur gewählt wurde, wodurch die Konzentration zu klein an- 
gegeben wird. In den vorliegenden Messungen unterscheiden sich 


- 
% 
Er 
ar 
1 
D 
5 ules: ‚= ıt dıeser Annahme läßt sıch dıe Konzentration der 
12 
F oe positiven lonen berechnen. Sie liegt bei den vorliegenden Messungen 
. > ständig um den Faktor 3+5 höher als die der Elektronen (bei sehr 
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f 
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die beiden Elektronengeschwindigkeiten um den Faktor 3; der Unter- 
schied zwischen n_ und n, würde somit beim Einsetzen der höheren 
Geschwindigkeit um Y3 größer, so daß wir schon ein Verhältnis 
von 5+9 für die lonen- und Elektronenkonzentration erhalten. 
Diese Zahlen sind durchaus in der Größenordnung, wie sie vielfach 
angegeben sind. Das noch vorhandene Überwiegen der positiven 
Ionen würde nach den neuesten Veröffentlichungen — insbesondere 
von Arnot — durch negative Ionen zu erklären sein, wie schon 
an anderer Stelle!5) ausführlicher angegeben wurde. 

In den Tabellen 7 und 8 sind für verschiedene Sondenlängen 5 
die Werte der Elektronen- und lonenkonzentration angegeben, Kr 
sowie das Verhältnis beider. 


Tabelle 7a und 7b 


d=08mm. n_,n,, = in Abhängigkeit von der Sondenlänge 
+ 


l n_ x10" n, x 10 n l n_ x 10% in, x 10" n_ 
mm em? em + mm cm? | em”? 
x | 
25| zı | 42,5 1,4 | 6,8 
5 | 116 44,0 3,6 11,3 
9 | 10,7 37,2 5,4 16,6 
13 10,0 36,0 45 17,2 
18,5 9,4 30,6 5,7 20,9 
24,5 6,5 25,6 6,5 25,6 
Tabelle 8a und 8 
=» Bas == in Abhängigkeit von der Sondenlänge 
- 
I jn_ x 10% in, x 10") n_ l n_ x 10" n, x 10%) n_ 
em3 em3 ny inm | em”? | | my 
5,7 64 lo 0,3 21 0,1 
9,7 67,1 (014 + +|i#35 1,6 91 [01 
15,9 73,9 /022 4 3,8 19,1 | 02 
18,5 841/022 5,8 298 | O01 
15,3 i* 40,6 | 0,2 
137 | 516 205) 90 | 511 | 
112 49,7 1023 .45| 11,2 49,7 | 02 
Fall 2 Fall 1 Be 
’ Durch Messungen mit einer verschiebbaren Sonde wurde der ie 
Einfluß des Sondenträgers bzw. der Gefäßwand auf den zur Sonde 


_ fließenden Strom untersucht. Bei langen Sonden zeigen sich in 


fii, 
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den halblogarithmischen Geraden Knicke, der Neigung des oberen 
Teiles entspricht die wirklich in der Entladung vorhandene Elek- 
-tronentemperatur. Bei kurzen Sonden wird diese Temperatur stets zu 
hoch gemessen infolge einer um den Sondenträger sich ausbildenden 
er Störschicht. Dieses ist darauf zurückzuführen, daß die Relaxations- 

längen für schnelle Elektronen von geringer Anzahl größer ist als 
für langsame Elektronen, daß also in die Störschicht eintretende 
‘schnelle Elektronen leichter zur Sonde gelangen können als lang- 
same und sich gegenseitig in der Störschicht weniger beeinflussen. 
Durch Messung mit kurzen Sonden an verschiedenen Querschnitt- 
stellen oder an derselben Stelle mit verschiedenen Bogenstromdichten 
: oy oder mit verschiedener Sondenlänge läßt sich auch die wahre Elek- 
_ tronentemperatur finden. Dieses ergibt sich aus der Boltzmann- 
5 schen Beziehung. 

Bei der Raumpotentialbestimmung gibt die Bestimmung des 
(oberen) Knickpunktes der In i_ — V -Geraden immer die richtigen 
Werte. Für die Rechnung als Zylinder- bzw. Kugelsonde muß die 
niedere wahre Elektronentemperatur eingesetzt werden. Die Ab- 
weichung bei den Bestimmungsmöglichkeiten liegt innerhalb der 
Meßgenauigkeit. Bei der Bestimmnng der Konzentration der Elek- 
tronen und positiven Ionen ergibt sich ein Überschuß der positiven 
Ionen um den Faktor 3—5, der hauptsächlich durch negative Ionen 
erklärt werden kann. 


rs od Die vorliegende Arbeit wurde in der Abteilung fiir technische 
sake und Akustik des Physikalischen Institutes der Tech- 

ischen Hochschule Braunschweig ausgefiihrt. Herrn Prof. Dr. E. 

Liibcke möchte ich fiir die Anregung zu dieser Arbeit und seine 
ständig fördernde Unterstützung und Herrn Prof. Dr. H. Diessel- 
horst für sein Interesse herzlichst danken. Zu den Messungen 
wurden Apparate benutzt, welche dankenswerterweise die Deutsche 
 Forschungsgemeinschaft Herrn Prof. Lübcke zur Verfügung gestellt 
hatte. 
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Eine vorläufige Mitteilung N 

über magnetische Nachwirkung 
Von Heinz Wittke 


(Mit 3 Figuren) 


Ich hatte in meinen beiden letzten Veröffentlichungen ') über 
‘oe Nachwirkung im Prinzip die in Fig. 1 skizzierte An- 
ee ordnung zugrunde gelegt. Darin bedeutete Sp eine Ringspule, die 


rhe 


) 


‘mit dem zu untersuchenden Ferromagnetikum gefüllt ist, ‘r ihren 
a _ Eigenwiderstand und R einen parallel geschalteten Widerstand, E, 
und E, die Spannungsquellen, von denen E, zur Zeit t= 0 durch 
den Schalter S abgeschaltet werden kann. Durch ein zur Zeit t an 
J ge die Sekundärwickelung von Sp gelegtes ballistisches Galvanometer 
“A u; wird die dann noch in diesem Moment in Sp herrschende Induktion 


b 
n=B+ mit b = 
bestimmt. 


Shin ist, wenn man von Hysteresis absieht. Dabei sprach ich die Ver- 
%: 3% mutung aus, daß das der Lösung (1) zugrunde gelegte Differentialgesetz 


(2) B=u-H+y(wH-— B) 


ARE 1) H. Wittke, Ann. d. Phys. [5] 28. S. 442. 1935, ausführlicher in Schrift. 
We, ws d. Kbg. Gelehrt. Ges. Naturw. Klasse 12. Jahr, Heft 3. 1935. x 


Mit dieser Anordnung konnte ich damals zeigen, daß für «> y 
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nur als Näherung zu betrachten sein dürfte und man wahrscheinlich 
nicht mit eenem y-Wert zu rechnen haben würde, sondern mit allen 


0=y= mw, wobei es Aufgabe weiterer Forschung sein müßte, durch ir ok 
Messungen über längere Zeiten t, als es mir damals möglich war ope 
(0 St < 1/,, sec), die Verteilungsfunktion der 7 zu ermitteln. 
0 0 Zeit 


4 


Va 
23 L 
46 


Fig. 2. 
Die Tab. 1 gehört zu den Kurven a, zu den Kurven b gehören diese Angaben: 
= 632,1 mOe; b/h = 439,7 Gauss/Oe; Ring Nr. 6v; Ringtemperatur +19° C; 


20. Januar 1936 

= Ich habe jetzt derartige Messungen für 0=t< 1,8 sec ab- 
geschlossen und dabei gefunden, daß Gl. ” zu ändern ist in BAe 

Betrachtet man also nur den Verlauf für größere Zeiten, so ist 


4 % t). 
(4 
Ein in Tab. 1 wiedergegebenes Meßprotokoll und seine und einer 


anderen Meßreihe graphische Darstellung in Fig. 2 sowie beider Aus- 
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Tabelle 1 


Am Anfange Am Ende Notierte GréBen 


der Messung der Messung 
if: ; 21,0 | 21,0 Galvanometertemperatur in + °C 
— 21,1 Zimmertemperatur . .in +°C | 
20,0 Ringtemperatur . . .in +°C 
2500 250,0 Eigenwiderstandr . . in Ohm 
im 250,0 | Parallelwiderstand R . in Ohm 
7% Mee 1,016 1,016 AuBenspannung E, . . in Volt 
2,02 2,02 Innenspannung E; . . in Volt 
y= 43,1 mOe. Shunt 3000. 15.8.1936. Ring Nr. 6n. 
Zeit in | Galvano- | 
are beliebigen | meteraus- | log b | n 
oa Einheiten | schläge | ™ fiir groBe Zeiten 
oo) 790 | 1.6987 7.070 ausschlag gemittelt aus 
ee 5,66 | 1,8455 8,352 5 Ablesungen bei + 7,6°/, 
DET, 9000 4,04 | 1,9899 9,885 größtemFehlereinerEinzel- 
ia 4 € 12000 3,16 | 2,0966 11,18 ablesung gegen das Mittel 
h m? Die Eichung ergab b/h = 354,3 Gauss/Oe. 


Behält man das Differentialgesetz (2) bei, so ergibt sich daraus 


€ 


unter Berücksichtigung von (3) als gesuchte Verteilungsfunktion der y 


fi) = at. yore. 

oe, Die mathematische Herleitung von (5) aus (2) und (3), die Messungs- 
ey ae ergebnisse über die Abhängigkeit der Konstanten ©, und a von b 
und von der Temperatur (und den Einfluß elastischer Verspannungen 


ct a auf sie) sowie die sich auf (2) und (5) griindende Diskussion der 
Zr 4 Messungen bei Wechselstromerregung statt der hier benutzten ape- 
riodischen bilden eine abgeschlossene Arbeit. Sie wird jedoch äußerer 
Gründe wegen erst in nächster Zeit in dieser Zeitschrift ausführlich 
abgedruckt werden. 


Königsberg Pr., I. Physikal. Institut, den 13. September 1008: 
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Über den magnetischen ‚Skineffekt 
in ferromagnetischen Kreiszylindern 
bei schwachen und starken Wechselfeldern’) 


Von Bruno Kyewski 


A. Einleitung 


Bringt man einen unendlich lang gedachten zylindrischen Leiter 
in ein longitudinales magnetisches Wechselfeld, so wird im Innern 
eine Wirbelstromverteilung auftreten, welche von dem Querschnitt 
des Leiters wie von der Frequenz und Amplitude des äußeren 
Wechselfeldes abhängt. Praktisch erzeugt man ein solches magne- 
| tisches Wechselfeld durch eine sehr lange Spule, welche von einem 
> Wechselstrom durchflossen wird. Gegenüber dem gewöhnlichen Skin- 
1 effekt haben elektrische und magnetische Feldstärke ihre Rollen ver- 
tauscht. Zur Berechnung der Feldverteilung im Leiter ist die Lösung 
der Differentialgleichung 


3 4H —j4npooH =0, 
w = Kreisfrequenz, 


o = spezifische Leitfähigkeit 
erforderlich. Nach der Integration liefert der komplexe Poyntingsche 
Vektor die Joulesche Wärmeentwicklung und die magnetische Feld- 
energie im Stab, dieselben sind für die induktive Heizung von 
Wichtigkeit. Für konstante magnetische Permeabilität läßt sich die 
obige Differentialgleichung mit Hilfe von Besselschen Funktionen 
lösen. O. Strutt?) hat dieses Problem eingehend behandelt. Hier 
soll die Rechnung für ferromagnetische Kreiszylinder durchgeführt 
werden, wobei man sich zunächst auf schwache Felder beschränkt. 
| Die Aufgabe wurde für den gewöhnlichen Skineffekt von E. Hinze *) 

näherungsweise gelöst. R. Becker‘) hat in einer Notiz zur Messung 
der Permeabilitat mit Hilfe des Hauteffektes einen Fehler der 
Hinzeschen Rechnung richtig gestellt. 
1) Göttinger Dissertation. 


2) O. Strutt, Ann. d. Phys. 82. S. 605. 1927; Arch. f. Elektrot. 19. 
8.424. 1928. 


4 R. Becker, Ann. d. Phys. (5) 27. S. 123. 


Hinze, Ann. d. Phys. (5]19. S. 143.1934 
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B. Magnetischer Skineffekt bei schwachen Wechselfeldern x 


Im Gebiet kleiner Feldstärken kann man nach Rayleigh die a 
Hysteresisschleife durch den Ansatz beschreiben: 


B = (u,+¢H,)H — za — für den aufsteigenden Ast, 
(1) 


für den absteigenden Ast. 


Beschränkt man sich auf die Frequenz der Grundwelle, so gelangt 
man durch te zur komplexen Darstellung: 


Die magnetische Induktion B und die Feldstärke H an sn 
um so stärker phasenverschoben, je größer die Amplitude ist. Nach 
Jordan besteht noch eine amplitudenunabhängige Nachwirkung; 
so daß man mit R. Becker den Ansatz zu: 


(8) 
vervollständigen kann. 


Die Differentialgleichung lautet für die magnetische Feld- 
verdrängung in dem Fall konstanter Permeabilität: 
d’H 1 = 
(4) oF yt. 


dr? r 


d’H 1 


\ 


Setzt man: 


so lautet die Differentialgleichung: BER 


dH 


Für die Lösung der Differentialgleichung sind zwei Fälle zu unter- 
scheiden: 
der Fall niedriger 


) der Fall hoher ‚Frequenzen 


3 
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a) Der Fall niedriger Frequenzen j= 
Zur Lösung der Differentialgleichung 3 
Ne 
(8) H=H,+H,+H,+---, 


wobei für alle Werte von r vorausgesetzt wird, daß |b| - | H 
Man erhält dann ein System von Differentialgleichungen: 


d’H 1 dH 


H,=0, 
d? H, 1 dH 
+ H, H,) H, 1) 
| 


Im Fall niedriger Frequenzen läßt sich die Lösung des Systems der 
Differentialgleichungen mit Hilfe eines Potenzreihenansatzes auf- 
finden, da 9 <1 und 0, < 1 sind. 


RR; 


Die nullte Näherung liefert dann’): 
H, =a-J,(¥-—G 0) 


(11) a 
36 - 64-64 


, 1-86 (48 — 46%) 


_ Die Entwicklung konvergiert gut für oe < 3, also in einem größeren 
_ Bereich als anfänglich verlangt. 

In Fig. 1 ist die magnetische Feldstärke H, = J,(V— (7 +8) - 9) 
als Vektor in der komplexen Ebene dargestellt, der Winkel g be- 
deutet die Phasendifferenz von H, gegenüber der Drahtachse. 


1) H ist konjugiert komplex zu H. 
2) Im = Funktion mter Ordnung. 
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e Da ? klein ist, berücksichtigt man bei der Näherungslösung nur 
Größen erster Ordnung in ß. 


(1+f*) (68+108%) 15 


(j +28) o?+ 


2* 


ST 


Vteil 

ER Fig. 1. Die magnetische Feldstärke in Abhängigkeit von 9 
j als Vektor in der komplexen Tun dargestellt 


= : ee Die erste Näherung in b ergibt: 

4° 36-64 ° 36-64-64 
es magnetische Feldstärke beträgt: 


a. Bre) (4-0 
2) 


3 
628 
4 
4 
| 
u te 
» 
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Für #=0 erhält man für die zweite ne in b: 


3) -H, =ba-Jal?- .02— 
(16) |H,| -H, =ba ‘Ts |; 
H, 2 105 74 

(18) 

ar... 


(3b-b) (9b-7h)j (6b+18D) 
2-64 9.36.64 36-64-64 
Die magnetische Feldstärke wird fiir 8 = 0: 

[ne — 32 36-64 36-64-68 

_ (80-4 Gb-7D5 (6b+18) 5 
| b-|a| (a ? 2.36-64 °° + 36.64.64 °° 


(19) H, = — ba-|a|*-| 


(20) 


In Fig. 2 ist H(Y—j-e) in Abhängigkeit von « dargestellt. Die 
Hysteresis bewirkt bei schwachen Wechselfeldern für «= 1 eine 
Vergrößerung von Phasenverschiebung und Amplitude der magneti- 
schen Feldstärke gegenüber dem Fall konstanter Permeabilitit. 

Der Konstante a bestimmt sich aus der Bedingung, daß die 
Tangentialkomponente der Feldstärke eine Grenzfläche zweier Medien 
stetig durchsetzt. 

H(o = o,) = H, = 0,4-n»n-t, 

wenn der Strom mit 7 und die Anzahl der Spulenwindungen pro 
Längeneinheit mit n bezeichnet werden. 


(21) J (V-GU + f)- 00) 
| b-| H,| % a = 64g. 36 +: 
IV = G+ B)+ + +B) 00), 


Der erzeugte Induktionsfluß a, 


Go fe 
0 


k? 
0 
oder 
( ) | | (28 16 Q*— ; 6-64 


4 
3 
| 
| 
4 
q | 
ha “inp 
| 
| 
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ir 9 =0 wird: 


| (eo: V-N-I V-)) 


20° _ 23° 690° 


16 6-64 8-36-64 10-36-64-64 


In Fig. 3 ist der InduktionsfiuB in Abhängigkeit von @ aufgetragen. 
Bei den Auswertungen wurden in den Reihenentwicklungen die 
Glieder bis zu o!° einschließlich berücksichtigt. 


0 OF 42 
a/teil 

“0,6 

0-10 
0-1 


Tmaginärteil 


Imaginärteil 


> 


Fig. 2. Die magnetische Feldstirke Fig. 3. Der Induktionsfluß 
als Vektor in der komplexen Ebene dargestellt 


Zusammenfassend läßt sich über die magnetische Feldstärke 
und für den Induktionsfluß (dazu vgl. man Fig. 8) sagen, daß die 
charakteristische Konstante « eine Vergrößerung von Amplitude und 
Phasenverschiebung bewirkt, 3 dagegen eine Erhöhung der Amplitude 
und eine Verminderung der Phasendifferenz. Die erzeugte Joule- 
(25) Q= 


an 


1936 
| 
3 
q | | 
4 | 0-05 | 
\ -0,5 
74 
N 
7 
56 7 \g! 
0# / A : 5 
Uw 
x 
Br Die erzeugte magnetische Feldenergie ergib 


‚schen Wärme und der 


B. Kyewski. Über den magnetischen Skineffekt usw. 
wobei F = Realteil (L) und G = — (L) ist und 
(a V y—)): J, (0% 


16 6-64 8-36-64 10-36-64-64 


Für « = 0 gehen die beiden Ausdrücke (25), (26) in jene Form über, 
welche M. O. Strutt in seiner Arbeit abgeleitet hat. In Fig. 4 
stellt F die dissipierte 42 
Wirme im Kern, G 
die Kernreaktanz als 
Funktion einer Wech- 
selstromkonstanten 
Kernhalbmesser dar. 
Die Kurven zeigen, 
daß in ferromagneti- 
schen Kreiszylindern 
(«> 0) bei schwachen 
Wechselfeldern eine ° 08 10 


Erhöhung der Joule- Fig. 4. Kurven für die magnetische Feld- Bait 
energie und Joulesche Wärme 


in Abhängigkeit von o dargestellt 


40 


magnetischen Feld- 
energie gegenüber dem Fall konstanter Permeabilität (2 = 0) eintritt. 


b) Der Fall hoher Frequenzen 
Da im Fall hoher Frequenzen der Skineffekt selır stark aus- 
geprägt ist, kann man die Oberfläche des Drahtes als eben be- 
trachten. (6) läßt sich auch schreiben: ar 


Für den aan hoher Frequenzen ist 9, >1. Führt man die 
Transformation s = 9, —@ ein, wobei man unter s die mit k mul- 
tiplizierte Entfernung von der Drahtoberfläche versteht (s < o,), S 
geht die Differentialgleichung über in: 

d? H ı 
de® ds 


(29) 

Beschränkt man sich auf Glieder nullter Ordnung in s/ Co» so fällt 
1 = 

8 

oberfläche als eine 


der Ausdruck — 


(30) —j(l—jf)-H H=0. 


7 


= 
631 
ae 
ie 
a — 
> 
4 
— 
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Man bestimmt eine Näherungslösung so, wie sie R. Becker durch- 
führte, indem man von b alle Potenzen höher als zweiter Ordnung 


vernachlässigt. 
Man setzt die Lösung in der Form an: 
(31) H=H,+H,+H,+--- 
H, soll klein von der Größenordnung b und H, von b? sein; 
dann wird: 


Da ß ebenfalls klein ist, so sollen bei der cn alle 
_ Quadrate und höheren Potenzen von $ vernachlässigt werden. 


Die nullte Näherung in b liefert: 


j-- 
(33) H,=C-e Ve ( 
Die erste Näherung in b ergibt: 


a Die zweite Näherung in b beträgt: 
e zwei a gi ) 


_jebee Vt 

: 2(2 +9) 

> 


+ j-0,331) + C3- 0,503] 2 


i Integrationskonstante C bestimmt sich aus der Bedingung, daB 
der Drahtoberfliche fiir s=0 H, = 04-+a-n-1% wird. 


= a : (0,508 + j - 0,331). 


1 


4 
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C=A.[ı 
a 3 5 4 
Tr 


Man erhält in erster Näherung für den erzeugten Induktionsfluß: 


2a 
0 


oo 


k? 
0 
1+j 


| Ve 4 , 


=|Cl-e— "(1 + v2 
4 ~ 

2ne(1 - — 

(41) O= | a tae, 


— Da .C.|e VE (1-5-3) 
(1+J-I1-j— 
2 


Va 
+ |C| @-(0,503 + j- 0,881) 
+ |C|?a?(0,184 + j-0,219) 2. 
e Ve a: 


+ (1 —j 1,006 + j-0,831) Y2- 


Die zeitliche Änderung des Induktionsflusses ist: 


(0,503 + j-0,381)a-|H,| + (0,182 + j-0,789)-«2-|H, |? +-- | 


(43) 


4 
u. In erster Näherung gilt für die Integrationskonstante C: EN 
. 
= 
245 
e 1 
a 
| 
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| 
ou 
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Die ra J Joulesche Wärme beträgt: 


(44) 


_ ER 503 em er 12 
(45) F = + 57. + 0,508 (2) |H + 0,182 2): 
Die erzeugte magnetische Feldenergie ergibt: 


(46) 


Al) = (1 — + )|H,| + 0,789 )?|H,|*+- 


Unter |H,| versteht man die Amplitude des Magnetfeldes an der 
Drahtoberfläche. 


70 


ao 


a 


Fig. 5. Die erzeugte ae Fig. 6. Die erzeugte 
Joulesche Wirme magnetische Feldenergie 


In Fig. 5 ist die dissipierte Wärme im Kern und in Fig. 6 
ist die Kernreaktanz fiir 9, = 2 bis 9, = 10 in Abhängigkeit von 


a = —- dargestellt. 
Ho 


me 


Bezeichnet man 
= Vu i+ „+ 0,508 A, + 0,182(*)". H |? + 


Vito 2u 9831 Fr HA, + 0,789 (4) - H, |, 
wobei man unter ur die wirksame Permeabilität der entwickelten 
Wärme und unter u, die der Selbstinduktion versteht, so läßt sich 
ein Näherungsausdruck linear in « für die wirksame Permeabilität u 


ableiten, wenn -|H,| <1 ist. 


ur = +7 + 1,006-¢- H, 
UL = My — + 0,662 - 5 
= 27 + -&-|H, 


zu a 
— 
. 
| 
5 
= 
RAS 
— ‘ 1:3 
2 
: A 
= 
% yn 
8 
WER 
T 
Er 
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Es ergibt sich also derselbe charakteristische Wert mit denselben 
Zahlenfaktoren wie bei R. Becker als Folge des Ansatzes fiir die 
der Permeabilitit. 


.C. Magnetischer Skineffekt bei starken Wechselfeldern 

_ Während der bisherige Ansatz für schwache Felder eine Er- 
nittiong der Größen y, «, somit der Kennwerte der Hysteresis- 
schleife, gestattet, soll jetzt der Fall starker Felder, welcher formal 
zum Teil in der bisherigen Rechnung erfaßt worden ist, nach dem 
summarischen Ansatz von Arkadiew berechnet werden. W. Arkadiew 
hat eine Theorie der Ferromagnetika bei hohen Frequenzen gegeben. 
Durch Einführung der magnetischen Leitfähigkeit o wird der durch 
die Phasenverschiebung bedingte Energieverlust berücksichtigt. 

Es wird gesetzt: 


(51) W“=u-—j- 
Der Phasenwinkel gw zwischen B und H bestimmt sich aus der 


Gleichung: 
(52) =i =). 
Die er B ist jetzt eine komplexe Größe, B wird ein Ausdruck 
von der Gestalt: 
(53) B=wWH=(u-j:o)-H. 
Im Gegensatz zur bisherigen Rechnung ist über die Feldabhängigkeit 
von g nichts angenommen, sondern man rechnet mit einem kon- 
stanten Mittelwert. Insofern ist diese Rechnung viel gröber. 
Mittelstraß!) hat rechnerisch und experimentell bei Untersuchungen 
über den elektrischen Skineffekt die Zweckmäßigkeit der Arka- 
diewschen Theorie geprüft. 
Die Besselsche Differentialgleichung läßt sich in diesem Fall 
streng integrieren. 
Die Differentialgleichung (4) lautet für komplexe Permeabilitat: 
(54) 


und dieselbe liefert die Lösung: 


Die A sich aus der Bedingung, daß die 
komponente der Feldstärke eine Grenzfläche zweier Medien stetig 
durchsetzt. Fürr=a wird: 


| 
e 
¥ 
4 
E 
3 
a a 
1) K. A. MittelstraB, Arch. f. Elektrotechn. 18. S. 595. 1927. a 
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| O4-a-n-isA- J, + b)) 

Der erzeugte Induktionsfluß ergibt: 


0 


_ - Vo G +d. 
| r-V (7 +) 


Die Berechnung der Jouleschen Wärme liefert: 


a-l-n?-@? 

| 

| F= Re 


Imoginärteil 


= Fig. 7. Die magnetische Feldstärke in der komplexen Ebene dargestellt 
fers: und für die magnetische Feldenergie erhält man: 
We 


wobei 
Die Berechnung der einzelnen Funktionen:') Die magnetische Feld- 
stärke ist durch J, (2) für i, sot 


z2=k-Y—(j+b) 


bestimmt. 


1) Vgl. Jahnke-Emde, S. 90/100. 


peu 


| 
= 
a 
m 
| 
— 
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4 Eine Reihe von Funktionen ist von MittelstraB berechnet 
worden. Diese Berechnung wurde noch durch einige weitere Kurven 
vervollständigt (vgl. Fig. 7). 

Der InduktionsfluB ist durch die Besselsche Funktion 
1. Ordnung bestimmt (vgl. Figg. 8 u. 9). 


(60) 


2 
für 2=k- V—(G +5) geht z- J, (2) ther 
@b-. (1 — 6b? + - 


(5b — 105° + b) . . [ 2bk* (1-35. ké 
(46 - 4 - ks 1- “ae 


Diese Reihe hat für k > 2,5 eine schwache Konvergenz, conte 
führt man semikonvergente Reihen ein. 


37 15 105, 
(62) V2:n»z j 82 128 z? 1024 2° 

5 


über in: 


e+ sin (0 + 


k-Y-(+b)-J, (k- Y- (j+b)= —_— 


- (cos 2 +7j sin f)- 


—g-sin(#+4) + 5+ 


- re 
(63) 105 sin 24 4725 cos (30 + 4) he ee 


| 15 cos (9 ++} 
—j- 00s + + 128-9 +- 


4725 sin (39+ 4 


ß = 0+c08 G + i) + 
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Fig. 8. Der Induktionsfluß in der komplexen Ebene dargestellt ö 
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für b= 0,0 wird 9 =k und der Ausdruck geht über in: 


4 


~~ 


2.J, (2) = 5 (cos 
4725 
+ 128kV 2 32 768k* 2 
15 105 4725 
128kV2 1024k* 32768k*/2 


Rea/tei/ 


zed, (2) 


um, so erhält man: 
J,(@) 


Rechnet man den Ausdruck 


— AT = — 0:sin + 


25 cos (so + 4 


. 


4 
3 
(64) | 
‘oe £ 
3 
23 
4 
> 
Fig. 9. Der Induktionsflu8 in der komplexen Ebene dargestellt 
3 Jo (2) J Be Be 12828 
Für E 
er f 
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Br Der Ausdruck für die Joulesche Wärme geht über in: 


(65) 


für b=0,0 und k=o 


k 
66 -" — 
(7 2 8kYV 2 128%? 2 


für sehr große Werte von k geht (47) über in: 


(67) (e - sin (# + 4) 


wird b=0,0 und o=k 


(68) 


Entsprechend erhält man für die magnetische Feldenergie: 


os (e+ cos (9 + 4) cos 24 
8.9 + 8. 0° 


(69) 


Für sehr große Werte von ” aa (51) über in: 


(70) = 008 (: +4)). 
Wird b=00 und o=k: 


W = 4 
Aus den Formeln (65) bis (71) folgt, daß die Ausdrücke für die J onlösche: 5 mae a 
Wärme wie für die magnetische Feldenergie für sehr große k-Werte 
in Geraden übergehen. Die spezialisierten Ausdrücke (68) und (71) 
sind jene, die M. O. Strutt in seiner Arbeit abgeleitet hat. 
In den Kurvendarstellungen (Figg. 10 und 11) sind die Joulesche 
Wärme und die magnetische Feldenergie für den Bereich k = 0,0 
bis k = 10,0 in Abhängigkeit von b= 0,0 bis b = 0,8 dargestellt. 
Vernachlässigt man bei der Entwicklung von a Y- ‘G +b) die 
höheren Potenzen von b: 


: fi 
B 
+ 
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5 
= 
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7 
% 
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0 


+? -b-j-Vyi + 


k 7 
= . — © 1 b? b, 
Q Vyi+ 


& 
= are 0 2 4 6 8 10k < 
Fig. 10. Kurven fiir die Joulesche Wiirme = 
| 
3 Fig. 11. Kurven für die magnetische Feldenergie 
| ER y 
| E 
4 


so erhält man eine lineare Beziehung in 0’: BE 
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(74) — 2-0’. 


berechnen kann. Nach einem Vergleich mit dem Ansatz von Ray- 
leigh-Jordan entspricht: 


D. Zusammenfassung 


Aus (51) folgt, daß man aus der Kenntnis von o’ und ü die Größe b 


1. Es wurden nach dem Ansatz von Rayleigh-Jordan für 
die Hysteresiskurve sowie nach dem von Arkadiew die Beziehungen 
für die Joulesche Wärme und magnetische Feldenergie berechnet 
und graphisch dargestellt. 

2. Die Funktion k-Y- (7 +b)-J, (k- Y—(j+b)) wurde von 
k = 0,0 bis k = 10,0 für b=0,1 bis b= 0,9 tabelliert und von 
k = 0,0 bis k = 6,0 graphisch dargestellt. 

3. Die Differenz von ur — u; ergibt nach Rayleigh-Jordan 
in erster Näherung eine lineare Beziehung in & und nach Arkadiew 
entsprechend eine für o’. 


Die Arbeit wurde im II. Physikalischen Institut der Universität 
Göttingen ausgeführt. Herrn Prof. Dr. Joos bin ich für die stetige 
Unterstützung und wertvolle Diskussion zu großem Dank verpflichtet. 
Herrn Prof. Dr. Reich danke ich für die Anregung zu dieser Arbeit 
wie für sein mir entgegengebrachtes Interesse. 


Göttingen, II. Physikalisches Institut und Institut für an- 
gewandte Elektrizität. 
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Der Zeemaneffekt der Alkalien 
nach der Diracgleichung') 


Von R. Schlatterer 

Mit Zeemans Entdeckung der magnetischen Beeinflußbarkeit 
der Spektrallinien im Jahre 1896 war die Atomphysik schon in 
ihren ersten Anfängen vor ein Problem gestellt, das ihre ganze 
weitere Entwicklung begleiten sollte. Zwar konnte die Lorentzsche 
Elektronentheorie sogleich eine Erklärung des Phänomens in seinen 
wesentlichsten Zügen geben, aber die späteren Beobachtungen über 
die Anomalien des Effekts und die an die Namen Paschen und 
Back geknüpften Erscheinungen in sehr starken Magnetfeldern 
gaben der Theorie immer neue Rätsel auf. 

Auch die Quantentheorie vermochte zunächst nur den „normalen“ 
Fall zu erklären. In der Deutung der anomalen Effekte kam ihr 
die Voigt-Sommerfeldsche Theorie zuvor, die im Grunde eine 
Vertiefung der ursprünglichen Lorentzschen Ideen war. Die Reich- 
weite der Voigtschen Vorstellung erstreckte sich auf die richtige 
Berechnung der Frequenzen wie auch der Intensitäten im Falle 
des D-Linientypus für beliebige Feldstärken. Uns sind so erstaunliche 
Erfolge heute kaum mehr verständlich, da wir doch an das Problem 
des anomalen Zeemaneffekts ohne die Kenntnis des Elektronenspins 
gar nicht mehr herantreten können. 

Die neuen Formen der Quantentheorie als Quanten- oder 
Wellenmechanik konnten ursprünglich von den Anomalien eben- 
sowenig Rechenschaft geben. Hierzu mußten immer erst irgend- 
welche Spinmodelle zusätzlich eingeführt werden. Solche Modelle 
lagen den quantenmechanischen Berechnungen des Zeemaneffekts 
von Heisenberg und Jordan?) wie auch den ersten wellen- 
mechanischen von Darwin?) zugrunde. 

Heute ist die von Dirac aufgefundene Wellengleichung des 
Elektrons der natürliche Ausgangspunkt für alle wellenmechanischen 
Berechnungen, in denen der Spin eine entscheidende Rolle spielt. 
Im Vergleich mit den spektroskopischen Tatsachen zeigt sich die 


1) Inauguraldissertation Universität München. re 
W. Heisenberg u. P. Jordan, Ztschr. f. Phys. 37. 8.270. 1926. 
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Überlegenheit der Diracgleichung gegenüber den älteren Ansätzen 
darin, daß sie die äußerst genau nachprüfbare Feinstrukturformel 
Sommerfelds in ihrer exakten Form liefert. Daß sie auch den 
Zeemaneffekt der Dublettsysteme in allen Einzelheiten richtig wieder- 
gibt, hat als erster Darwin’) gezeigt. 

Die vorliegende Arbeit will ebenfalls eine Ableitung der Fein- 
struktur und des Zeemaneffekts aus der Diracgleichung geben, und 
zwar für deren genauen Gültigkeitsbereich: bei Atomen mit einem 
äußeren Elektron. Als Ausgangspunkt ist dabei nicht die lineare, 
sondern die iterierte Diracgleichung gewählt, von der Sommerfeld 
und Maue?) gezeigt haben, daß sich ihre Lösung im Falle eines 
beliebigen stationären elektrischen Feldes an die der Schrödinger- 
gleichung anschließen läßt. Die Methode ist ein Näherungsverfahren 
und stellt sich als eine Entwicklung sowohl der Eigenfunktion als 
auch der Eigenwerte nach Potenzen der Feinstrukturkonstanten « 
dar. Dieser Näherungsansatz darf wohl für die Behandlung der 
Spektren als der der Natur der Sache am besten angemessene be- 
trachtet werden. Gegenüber der Paulischen Näherung, die eine 


Entwicklung für kleine Werte von 3 = = ist, ist der Unterschied 


analytisch durch das Auftreten der Biquaternionen an Stelle des 
Quaternionenkörpers gegeben. 

Im Gegensatz zu den bisherigen Berechnungen) des Zeeman- 
effekts nach Dirac, die die magnetische Störung an die exakte 
Lösung der Diracgleichung ohne Magnetfeld anschließen, werden 
wir die Feinstruktur- und die magnetische Aufspaltung gleichzeitig 
als Störungen behandeln. Dies ist konsequent, da ja die beiden 
Effekte tatsächlich von derselben Größenordnung sind. Es wird sich 
zeigen, daß in unserem Falle das Vorhandensein eines a 
kaum eine Komplikation der Rechnungen bedingt. I 2 


: + § 1. Erläuterung des Näherungsansatzes und Integration 
7 der ersten Näherung 


Das Leuchtelektron bewegt sich in einem zentralsymmetrischen 
elektrischen Feld mit dem Potential V = V(r), dem ein homogenes 
Magnetfeld 5 = §, mit dem Vektorpotential 

(1,1) ye 


1) C. G. Darwin, Proc. Roy. Soc. A 118. S. 676. 1928. 
2) A. Sommerfeld u. A. W. Maue, Ann. d. Phys. [5] 22. S. 629. 1935. 
(Im folgenden mit a. a. O. zitiert.) 
3) Vgl. außer C.G. Darwin (a. a. 0.) etwa L. de Broglie ,,L’électron 
S. Sauter, Stocks. f. 63. 811. 1930. 
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2) 


ee (mit etwas veränderter Bedeutung des Energieparameters, vgl. a. a. O. 


Eigenfunktionen für ein elektrisches Zentralfeld. 


4 Größenordnung «* zuschreiben, ist für den Feldstärkenbereich des 
' anomalen Zeemaneffekts deshalb gerechtfertigt, weil hier die ma- 
__ gnetische Energiestérung von der Größe der Feinstrukturaufspaltung 


ea _ Naherung von der Ordnung a? ist. Letzteres ist im Falle des 
_ Wasserstoffs von vornherein aus der Sommerfeldschen Fein- N 
strukturformel zu entnehmen und der Größenordnung nach gilt eS, 
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überlagert ist. Wir gehen aus von der iterierten Diracgleichung 
des Problems 


— Vy — eW— E,2w 


1 
dw+ 


| = qe (7 grad V) y+ 25° grad y) + 5° w. 

Mit y und yy bezeichnen wir die Vektoren 


> > 
wobei noch 7,7, = 7,, gesetzt ist. % : 
In der Gleichung (1,2) sehen wir die mit 7 behafteten „Spin- 


glieder“ und ebenso die vom Magnetfeld abhängigen Terme als kleine 
Störungen an. Die ungestörte Gleichung ist eine Schrödingergleichung 


8. 630) und hat als Lösung die bekannten Schrödingerschen 


Um aus der Schrödingerschen Lösung als nullter Näherung 
eine angenäherte Lösung von (1,2) zu gewinnen, setzen wir für die 
Kigenfunktionen yw und ebenso die Kigenwerte E’ Reihenent- 

 wieklungen an 


(1,3) 


die nach Potenzen der Feinstrukturkonstanten @ fortschreiten. Es 
soll also w, von der Größenordnung a" yw, und &, von der Ordnung 
a" FE sein. In der Gleichung (1,2) selbst betrachten wir als klein von 


1. Ordnung in & das Glied: 7 (7 gradV), 


2. Ordnung die Glieder: 3 (A grad); ie (758). 


Daß wir den vom Magnetfeld linear abhängigen Gliedern die 
ist, welche ihrerseits gegenüber der Rydbergenergie als nullter ia : 


_ das gleiche auch für die Alkalien. Das Glied e? %?/A*c? wäre rae. 
_ konsequenterweise als klein von vierter Ordnung anzusehen und . 
oe: bleibt daher fiir uns ganz auBer Betracht. 


2 
| 
! 
Bi 
ur 


Aus den Ansätzen (1,3) erhält man durch Einsetzen in (1,2) 
und Zusammenfassung aller Glieder von gleicher Größenordnung die 
Gleichungen, denen die aufeinanderfolgenden Näherungen genügen 
müssen. Für , ist es die (homogene) Schrödingergleichung: 


(1,4) Dy = (E*— 2EV— = 0 
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Ihre Eigenwerte und Kigenlésungen dürfen wir als bekannt 
annehmen. Im allgemeinen Fall, wenn E entartet und ein p-facher 
Eigenwert ist, gibt es dazu p linear unabhängige Eigenfunktionen, 
die wir durch einen oberen Index u unterscheiden: 

(1,5) Dwe=0 für u=1,2,---,P. 
Wir wollen voraussetzen, daß diese yw gemäß der Gleichung 


(1,6) [re V) ped E, 


orthogonalisiert und normiert seien. ee möge immer den 
Vorzeichenwechsel der imaginären Einheit andeuten. Die hier ge- 
wählte Form der Orthogonalitätsbedingung ist daraus zu verstehen, 
daß der Eigenwertparameter E in (1,4) quadratisch vorkommt (vgl. 
a hierzu die analoge Orthog. Bed. a. a. O. Gl. (14) für = E). 
Die allgemeinste zum Eigenwert E gehörige Lösung von (1,4) 
ist gegeben durch eine Linearkombination 


(1,7 = a. wit 
) Wo 2% w 


Die erste Näherung w, unseres Problems muß 
 Schrödingergleichung 


(1,8) Dw = grad V) — c? 


erfüllen. Die Konstanten «, und a, sind _ den allgemeinen 
Vorschriften der Schrédingerschen’ Störungstheorie aus der Lös- 


(1,9) f w"Dwdr=0 für 


zu Es muß also gelten: 


(7 grad V) — 2 


| D'a, 


- 
As >46 
“) 
¢ 
4 
Bu 
. 
Bas 4 2 
4 
— 
fw ( y grad V wedra,. 
mora: 
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Wir wollen zeigen, daß das restliche Integral gleich Null ist?). 
Nach den Gl. (1,4) und (1,5) ist nämlich te n Pe 

0 = grad (D we) — grad yr D yw 


= wm A grad — grad yw 4 ym — m ve. 


Integrieren wir diese Gleichung über den unendlichen Raum, — 


so ergeben die beiden ersten Glieder nach dem Greenschen Satz ER 
ein Oberflächenintegral, das aber wegen des raschen Abfalls der cS Ces 
Eigenfunktionen nach außen verschwindet. Wir erhalten daher: er iS 
(1,11) 0 = [ grad 

Damit ergibt sich aber aus (1,10): 

O=——" 4,4, fir m =1, 2,---,p. a 

Es folgt, daß notwendig «,= 0 sein muß, während die a, un- 

bestimmt bleiben. Die erste Näherung reduziert sich somit auf die I 
Gleichung: 
(1,12) Dy, = —y (7 grad 2) a we. 
Thre allgemeine Lösung ist a. a. O. [GL (22)] a Wir aS 
schreiben sie in der Form): 


(1,13) 12 (7grad I) a, w) + 2 b,wr. 
u 


1) Wie a.a.O. S. 633, Anm. 1 erläutert ist, haben wir dies als eine een ae 
wellenmechanische Formulierung des Impulssatzes aufzufassen. In der dort ean 
angegebenen Gi. (20) ist wieder « = E zu setzen. ER 
2) Man kann sich in folgender Weise sehr einfach davon überzeugen, ES x 


daß (1,13) das allgemeine Integral von (1,12) ist. Aus Gl. (1,4) schließt man: 


2E 
grad Dy, = D grad yw — 


> 
Durch skalare Multiplikation mit dem Vektor sen y folgt weiter: 


Im Verein mit (1,7) zeigt diese Relation, daß 
5 
ein partikuläres Integral der Gl. (1,12) ist. Ihr allgemeines 
Integral entsteht durch Hinzufügen einer beliebigen Lösung der homogenen ee 
Gl. (1,4). 


Ich verdanke diesen Hinweis Herrn Dr. Welker. 


ch 
| 
4 
a 
| 
| 
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Dies stellt bei beliebiger Wahl der Koeffizienten b, eine Lö- 
sung der Gl. (1,12) dar. Auf die nähere Bestimmung der b, brauchen 
wir nicht einzugehen, da sie in den folgenden Rechnungen keine 


Rolle spielen werden. 
; § 2. Das Störungsschema für die zweite Näherung u 


Da die Energieniveaus in erster Näherung nicht gestört werden, 
nennen wir jetzt die Energiestörung zweiter Ordnung einfach & 
statt &. Sie ist zu berechnen aus der 2. Näherung der Gl. (1,2), 
nämlich: 


| + A: + he 


In das erste Glied der rechten Seite setzen wir fiir w, die 
Darstellung (1,13) ein und führen in dem so entstehenden Ausdruck 
die folgende Umformung durch: 


own (grad grad > a, we) _ (r7 [grad V, grad >’ a, wel) 

7 u 
Daß hierbei die Größe y, herausgefallen ist, wird uns für die weitere 
Rechnung von Vorteil sein. Wir können also jetzt Dy, ganz durch 
die Grundlösungen w* ausdrücken: 


Da (grad V, grad - 


( y grad J ) > bw + 2 — | grad > w) 
+ | Da, w. 


Die Grundlage fiir die Stérungsrechnung bildet wieder die Lés- 


Was 
“7 
— 
7» 
: 
4 
a 
= 
= 
| 
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welche für die Werte m = 1,2, ---, p erfüllt sein muß. Setzen wir 
hier die vorangehende Darstellung für Dw, ein, so zerfällt (2,3) in 
eine Summe von Integralen. Einige davon ne wir ausführen: 


(2,4) > a, we dr = um» 


“ 


Zufolge der (1,11) ist ferner: 
, (2,5) fw un (Farad y) 2» wdr=0. 


Da das Magnetfeld 9 konstant int, wird auch der Wert des _ 
Integrals 


[FEGO 2G) Da 
durch die Gl. (1,6) geliefert. Es ist nämlich bis auf Glieder von der = 


da sowohl E—E, als auch V von der Ordnung der Belbuupingie 
(= «®E,) sind. Daher wird mit einer korrekten Vernachlässigung 


dun= [ E-V)wdr=E,[ 

und infolgedessen 

26) Dia, war = 779) Sadun. 


Wir führen nun noch die Abkürzungen ein: 

um + Vi wdr—2ie 

(2,7) ) 2 
| (rx, (grad V, grad 

(2,8) Run (Se) V, grad 


R,.m ist ein räumlicher Vektor, dessen rechtwinklige Kompo- 
enenten X,m, Yum, Zum Wir uns den entsprechenden Komponenten 
des unter dem Integralzeichen stehenden vektoriellen Produkts zu- __ 
geordnet denken. 


> 


a 


| 
> 
Agrad we) dr 
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Für die Auflösung dieses Gleichungssystems ist es offenbar 
nötig, hyperkomplexe Werte der Unbekannten a, zuzulassen. Natür- 
lich ist es wünschenswert, von vornherein einen möglichst einfachen 
Ansatz zugrunde zu legen. Wir bemerken daher, daß in (2,9) 
explizit nur Zahlen der Form 
(2,10) 7,9) 
auftreten, wo q und q komplexe Größen sind. Diese Zahlen bilden 
einen Körper (nämlich den der Quaternionengruppe). Es genügt 
deshalb sicher, die a, ebenfalls als Zahlen dieses Körpers anzu- 
setzen. Man wird sich aber weiter den Umstand zunutze machen, 

daß der Zahlkörper (2,10) noch Nullteiler!) enthält, mit deren Hilfe 
man durch rechtsseitige Multiplikation die Zahl der unabhängigen 
Parameter von vier auf zwei reduzieren kann. Ein solcher Null- 
 teiler ist zB. 1+%y,,. Wir machen also für a, den Ansatz: 
145) + *%23(— + 
= (Z, + 4733 re! + ay 13)° 
Man gewinnt hierzu leicht die Multiplikationsregeln 
Ya: 4, = (— vy, + + | 
19, = Y, — + 172)» 
=(—iz, — Y,)(1 

Indem wir tng Ansatz fiir a, in (2,9) eintragen und die vor- 

stehenden re benützen, bekommen wir die Gleichungen: 


Dieses sich mit komplexen Werten 
von z,, y, erfüllen, wenn wir setzen: 


> |(B..- i2.. + = 0 


“ 


h 


u 


(2,12) 


1) Vgl. hierzu etwa W. Franz, „Zur Methodik der TG 
Sitz.-Ber. d. Bayr. Akad. d. Wiss. (1935), 


50 
| 
’ 
3 
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Damit haben wir für die 2p Unbekannten z , y, ein System WE: 
von 2p linearen homogenen Gleichungen gefunden. Die Lösbarkeits- u : 
bedingung des Systems, nämlich das Verschwinden seiner Deter- N 
minante, wird uns in bekannter Weise die Energiestörung « liefern. Er 

§ 3. Berechnung der Termaufspaltungen RE 
Für die wirkliche Berechnung der Energi>störungen & aus dem ap “a5 


Gleichungssystem (2,12) ist es nötig, auf die darin vorkommenden 
Größen Bia, Zum etwas näher einzugehen. 
Die Eigenfunktionen nullter Näherung sind nach Schrödinger 


(in Polarkoordinaten) 


1 
(3,1) yt = R,,, (r) P,* (cos 9) ei“r. 


Sie erfüllen die Orthogonalitätsbedingung (1,6) nicht exakt, 
aber wegen der schon benützten Beziehung E — V = E, mit der 
hier erforderlichen Annäherung. Als Definition der Kugelfunktionen 
wollen wir die von Darwin eingeführte 

d \t+u (uy? 


zugrunde legen. Sie führt auf das Normierungsintegral 


tn 


(3,3) f [P, (cos sin = = 
0 
Außerdem nehmen wir an, daß 
(3,4) (r)rdr = 


ist, so daB also auch die Normierung der Funktionen (3,1) mit (1,6) 
übereinstimmt. 

Die Entartung liegt in den 2!+1 der 
magnetischen Quantenzahl 


(3,5) p=—l, ---, +1, 


die zu einem bestimmten Energieniveau E = E,, gehören. Wir 
denken dabei speziell an die Verhältnisse zu den Alkaliatomen. 
Der Fall des Wasserstoffatoms, wo sich die Entartung auch auf 
die azimutale Quantenzahl / ausdehnt, würde eigentlich im folgenden 
etwas umständlichere Überlegungen nötig machen. Deren Durch- 
führung hätte aber keinerlei Besonderheiten zur Folge. Im Er- 
gebnis haben wir vom Standpunkt unserer Rechnung das Wasser- 
stoffatom als ein Atom des Alkalitypus aufzufassen, bei dem ledig- 
lich gewisse Energieniveaus zusammenfallen. Auch im Magnetfelde 


Fs 
hs 
FN, 
> 
U. 
| 
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spielt die Entartung in / keine Rolle, während die u-Entartung von 
entscheidendem Einfluß ist. 

In den in (2,7) und (2,8) eingeführten Größen lassen sich nun 
die Integrationen wenigstens für die Winkelanteile alle ausführen. 


Da die Orthogonalität der “ von der besonderen Form der Ra- 
dialteile R,, nicht abhängt, erhalten wir: 


=A 


grad yr) de =f wy inHyndr = mHö,,.- 
dV @ 

4 (grad V, grad = [ ws SE dr 


dr + (Fe) dr dr 
SER 


Das in dem ane (2,8) vorkommende Vektorprodukt hat in 
Polarkoordinaten die Komponenten 


dV 1 éy* 
[grad V, grad we], = 0, [,=— sane 


(I, = _4V 1 dg 
dr r 


1 dV dy* 
(3,12) [].=[],cos # —[],sin = — —— 


Setzt man dies in (2,8) ein, so ergibt sich fir die reer 
nente von &,,,, 


2 
4 
: 
3 9 B A h c H 2 
( ’ ) um e+ +™ 2E, um? 
; 2000 wobei A nur noch von den Quantenzahlen n und / abhängt: hr 
nl 2 
R,.r dr. 
Pes a. ür die kartesischen Komponenten erhält man daraus: 
Mgt, 
(3,11 
| 1 dV 
zZ [ru dr = ( ) | yw" — —— iu wedr 
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wobei 


gesetzt ist. Entsprechend berechnen wir für die z- und y- -Kom- 
ponente 


SP [= ucotg  P® + 


u 
der Definitionsgleichung (3,2) gewinnt man durch Difie- 
(3,16) — pootg  P* + —! =— Pr 


mit deren Hilfe man auch das in (3,15) noch verbliebene Integral 
über # auswerten kann: 


er Nm Nip 


0 

2ni 1 2 

+ mt — m)! 


Das Ergebnis ist also: 


(It: 


Setzen wir noch zur Abkürzung 

(3,20) a= — iQ Wu: 
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<4 
be 

EN 
(3,15) 
| 
[| 
| 
x 
L \? 112 
} \ 
1 \2 ER 
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Die Beziehungen (9), (13) und (20) dieses Paragraphen führen 
wir jetzt in das Gleichungssystem (2,12) ein. Nach Ausführung der 
erhalten wir: 


le + A +mQ + (m— 1) m— Wa - 30 = 0, 


(3,21) 


— + [e+ A—mQ + (m+ eH|y,=0. 


Diese Gleichungen miissen fiir alle ganzzahligen Werte m von 
 — Ibis +1 erfüllt sein. Um die Eigenwerte « zu finden, für 

welche das System Lösungen besitzt, betrachten wir zunächst die 
a Spezialfälle m = +1. Da nach (3,19) 

W_.-ı1 = Wi, =0 

ist, haben wir die Gleichungen: 


TE 


ER In den übrigen Gleichungen des Systems (8,21) kommen die 
Be - nbekannten x_, und %, nicht mehr vor. Daher gewinnen wir aus 


E den Gl. (3,22) und (3,23) die beiden Eigenwerte 

er ite Sie hängen linear von der Feldstärke H ab und liefern fiir 
den s-Term (J = 0) die Aufspaltung 


+ 


Wir stoßen hier auf „Spin: wegen 
deren die einfache normale Zeemanaufspaltung nur bei den Sin- 
gulettsystemen in Erscheinung treten kann. ie 
ane In den übrigen Gleichungen des Systems (3,21) ist jetzt 
!+1=m=l der Wertebereich von m für die erste Zeile und 
der Wertebereich für die zweite Zeile. Wir be- 

 seitigen diese Unsymmetrie und werden zugleich auf den einfachsten 
Weg zur Eigenwertbestimmung geführt, wenn wir uns die spektro- 
_ skopische Erfahrung zunutze machen und an Stelle des ganzzah- __ 
ligen m eine  halbzahlige einführen. Dazu ersetzen wir 


& 
= 
| 
| 
Dies ist das Doppelte der normalen Zeemanaufspaltung 
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in der ersten Zeile von (3,21) den Index m durch m +}, in der 
zweiten m durch m — 4. Danach hat man 


1 1 

| [e+ A + (m + 3)@+(m- 5) Tm +h — = 0, 

| + +4-(m-5)0+ 
als neues Gleichungssystem, wobei nun in beiden Zeilen der Index m 
alle halbzahligen Werte von — ( _ 3) bis + (1 _ 5) durchläuft. 
Man bemerkt, daß die aus (3,27) für verschiedene Werte von m 
sich ergebenden Paare von Gleichungen keine Unbekannten gemein- 
sam haben. Daher erhält man auch für jedes m eine besondere 
Lösbarkeitsbedingung, nämlich die Determinantengleichung: 

1 1 


Die Lösungen dieser quadratischen Gleichung für « sind: 


829) tor 


Dies gilt zunächst für den Wertebereich |m| =1 — = Man sieht a 
aber nachträglich, daß diese Formel auch die Spezialfälle (3,24) RAR, a 
umfaßt, wenn man für das positive Vorzeichen der Quadratwurzel Sa i 
noch die Werte m = + ( + 3) zuläßt. Wir fassen daher das Er- Ba, : 
gebnis unserer Rechnungen für die Eigenwerte der Energie zu- LEN 


sammen in die Formeln: 


für 


2 
Die Feinstruktur der Spektralterme von Atomen mit einem 
äußeren Elektron wie auch ihre Veränderungen in starken und 
schwachen magnetischen Feldern sind hierin enthalten (§ 4). Ferner 
können wir auf Grund der Gl. (3,27) zu jedem Eigenwert die Koeffi- 
 gienten z, y, und damit nach (2,10) und (1,7) die zugehörigen 
Eigenfunktionen y, ermitteln. Deren Kenntnis gestattet uns weiter- 
hin die Berechnung der Auswahlregeln und Intensitäten der Spektral- 
linie 


> 
A 
| 
‘a wit 
| 
1 1 2 1 
| é=—A— mo — 3) Q?—mQea+ 
2 
| 
| 
| 
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§ 4. Diskussion der bisherigen Ergebnisse 


Die Zweiteilung der Gl. (3,30) bringt die Dublettstruktur der Alkalitern:; 
zum Ausdruck. Die erste Zeile liefert im Magnetfeld 21 +2 Werte für s. 
= statistisches Gewicht ist also 21+ 2, was einer inneren Quantenzahl 


=1+ 1 entspricht. Die zweite Zeile ergibt nur 21 verschiedene Werte 
ae für e, entsprechend einer inneren Quantenzahl j =! — -. 
Wenn das äußere Feld verschwindet, haben wir in (3,30) @ gleich Null 
zu setzen. Es wird dann: 


&=—-A+/1Q Sir Tom 1+, 


2 


(4,1) 
2 


Dabei ist aber zu beachten, daß wegen der in (3,30) angegebenen Ein- 
schränkungen über den Wertebereich von m die zweite dieser Gleichungen 
_ wegzulassen ist, wenn / gleich Null ist. Der s-Term ist daher immer einfach. 
Als Dublettabstand hat man nach (4,1): 


As = (21+1)Q. 


Die Bedeutung der in (3,14) eingeführten Größe @ ist hierdurch ge- 
geben. Ihre Berechnung würde die Kenntnis des inneren Atomfeldes V= Vir) 
und der zugehörigen Radialteile R,„. (r) der Eigenfunktionen zur Voraussetzung 
haben. Um jedoch die aus der Spektroskopie bekannten allgemeinen Regel- 
_ mäßigkeiten') zu verstehen, die in der relativen Lage der Dubletts und in der 
Abhängigkeit ihrer Abstände von den Quantenzahlen n und I hervortreten, 
haben wir eine solche exakte Kenntnis nicht nötig. Dies ist vielmehr aus 
einer qualitativen Betrachtung der Darstellung 
_ he fa V 

unschwer zu entnehmen. 

Wir sehen zunächst, daß die Größe Q immer positiv sein wird, da ja 


das Leuchtelektron überall unter dem Einfluß eines Anziehungspotentials steht 
und daher dV/dr für alle Werte von r größer als Null sein muß. Nach (4,1) 


%-4° 


Der statistisch bev orzugte =) Term liegt demnach auch energetisch 


durchweg höher als der (x - 9 3) Term. Da das statistische Gewicht der Terme 


für die Intensitäten der zugehörigen Linien und ihre Energie für die Fre- 
quenzen der Linien bestimmend sind, findet die Regel (4,4) zwangsläufig ihren 
Ausdruck in dem regelmäßigen Intensitätswechsel der Dublettserien der Alkalien. 


1) Vgl. A. Sommerfeld, „Atombau und Spektrallinien“ I, Kap. 7. 


1 =a 
&=-A-(l+1Q » » J=l=3. 


656 
_ 
= 
v 
(4,2) 
A 
4 
186 also stets 
= 
Pr 
= 
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Uber den Verlauf des vom Atomrumpf herrührenden Potentials V (r) 
können wir wenigstens folgende Aussagen mit Sicherheit machen: In der Nähe 
des Kerns ist dessen Anziehung allein maßgebend und daher 

Ze: 
r 
am Kern. Außerhalb der Ladungswolke ist die Kernladung durch die Rumpf- 
elektronen bis auf ein Elementarquantum abgeschirmt, also a: 


2 
r 
an der Peripherie. Wir fassen beides zusammen in die Formel: 
2 


wobei { eine nicht näher bekannte Funktion von r ist, die sich aber fir r > 0 
dem Wert Z und für r> © dem Wert 1 nähert und die für mittlere r 
jedenfalls innerhalb dieser Grenzen bleibt. An Stelle der Gl. (4,2) können 
wir somit schreiben: 


co 
2 WE 
; 


Man sieht, daB wegen des Faktors 1/r im Integranden in erster Linie 
der Verlauf der radialen Eigenfunktionen in Kernndhe fiir den Wert dieses 
Integrals maßgebend ist. Da nun die Aufenthaltswahrscheinlichkeit des 
Leuchtelektrons in der Nähe des Kerns um so geringer wird, je größer n 
ist, ist es klar, daß mit zunehmender Hauptquantenzahl die Dublettabstände 
kleiner werden miissen. Die durch das Schema 


NPs 
7 


(4,7) 
n 
Py 


np, — np, = Dublettabstand 
? 


charakterisierten Hauptserien der Alkalien spiegeln diese Regel deutlich wieder. 

Auch die Abhängigkeit der R,; von der azimutalen Quantenzahl / ist 
uns in großen Zügen bekannt. Den kleinen Werten von / entsprechen die 
„Tauchbahnen“, d.h. also große Aufenthaltswahrscheinlichkeit in Kernnähe 
und somit nach unseren Überlegungen großer Dublettabstand. In den Zu- 
ständen mit größeren ! hingegen bleibt das Elektron an der Peripherie des 
Atoms: also hat das Integral (4,6) einen kleinen Wert und damit auch der 
Dublettabstand. Spektroskopisch zeigt sich dies darin, daß in der Regel schon 
in der I. Nebenserie mit dem Linienschema 


sp, nd, 
3 ? 


die d,-Differenzen (l = 2) unmerklich sind gegenüber den p,-Differenzen (/ = 1). 

Wir wollen nun die Aufspaltungen der Alkaliterme im Magnetfeld auf 
BER Grund der Formeln (3,30) untersuchen. Entsprechend der Beziehung (4,2) 
: ie werden wir zunächst die darin vorkommende Größe Q durch den Dublett- 
abstand ersetzen und das Verhältnis 
de, _ Dublettabstand 


® normale Zeemanaufspaltung 


4 REN. 


| 
| 
| 
| 
| 
| 
| 
| 
| 
FL 
| 
3 
| 
| 
i 
in 
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einführen. Wir erhalten: 


1 
1 ‘ + fir jel+s, 
’ 


1 
— fir 


Formel') gibt den der Dublettsysteme mit Einschluß 


ne mit der Erfahrung. Sie wurde für den etwas spezielleren Fall 
des D-Linientypus in dieser geschlossenen Form erstmals von Sommerfeld?) 
im Anschluß an die Voigtsche Theorie des Zeemaneffekts abgeleitet. 

k Die Abhängigkeit der Termgröße von der magnetischen Feldstärke ist 
danach im allgemeinen keine lineare. Vielmehr sind die durch (4,10) gegebenen 
&, H-Kurven Hyperbeln. Die Abweichungen vom linearen Verlauf sind jedoch 
7 bei den meisten Termen innerhalb beträchtlicher Feldstärkebereiche ziemlich 
geringfügig. Man kann daher die exakten Formeln vielfach durch übersicht- 


1 + fir j=1+5: 
(4,11) .-4=-[1: - ) mo 
- fürj=1-—- 
2 
1 
Diese Beziehung ergibt sich als Spezialfall mit S = > L=1,J=1+ a 


aus der allgemeinen Landéschen g-Formel: 
J(J 
(J + 1) 
welche zwar auch nur für schwache Felder, dafür aber für beliebige Multi- 
plizitäten gilt. 

Ist andererseits das Feld so stark, daß die Feinstrukturabstände klein 


werden gegenüber der magnetischen huligaiteng, so werden \ wir in (4,10) @ 
und v vernachlässigen: 


I 
1 + fir j=1+ +> 
(4,12) j= -m+5)o 

— fir j=l- > 
2° 
Wir bemerken, daß die für das untere Vorzeichen entstehenden, vom 


1 - 
( - 3) -Niveau herrührenden, Terme aufgehen in den vom Niveau j = | + 
stammenden. Diese starke Vereinfachung des Aufspaltungsbildes bei as 
Feldstärken ist für den Paschen-Backeffekt charakteristisch. 
E Den vereinfachten Formeln (4,11) und (4,12) für die beiden Grenzfälle 
ist gemeinsam, daß sie die Größe Q nicht mehr enthalten und daher auch 


Bar von der Hauptquantenzahl n unabhängig sind. Dies stimmt überein mit der 


1) Sie stimmt überein mit der in „Atombau und Spektrallinien“ I S. 532 
angegebenen. Das dortige M ist unser — m, was lediglich eine andere Nume- 
 rierung der Terme bedeutet. 

2) „Zur Voigtschen Theorie des Zeemaneffektes“, Gött. Nachr. 1914, 


wi 
=» 
. 7 RS lichere lineare Entwicklungen ersetzen. So erhält man aus (4,10) im Grenzfalle 
2 schwacher Felder durch binomische Entwicklung: 
waa 
yx! 
. 
a 
= is 
3 
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empirisch gefundenen Prestonschen Regel'), wonach alle Linien einer Serie 
den gleichen Zeemaneffekt zeigen. Hiermit ist zugleich der Gültigkeitsbereich 
dieser Regel abgegrenzt. 

Die Feinstruktur des Wasserstoffspektrums ist vom Standpunkt der 
Diracschen Theorie schon so oft behandelt worden, daß wir hier nur kurz 
anzudeuten brauchen, inwiefern dieser Fall in unseren Rechnungen enthalten 
ist. Da für ein Coulombfeld auch die Radialteile R,, der Eigenfunktionen yw“ 
bekannt sind, können wir hier die Größen A und Q nach den Gl. (3,10) und 
(3,14) explizit berechnen. Es zeigt sich, daß beim Fehlen eines äußeren Feldes 
die besondere Form dieser Größen bewirkt, daB gewisse bei den Alkalien ge- 
trennte Energieniveaus zusammenfallen. An Stelle der Gl. (4,1) haben wir hier: 


Die Terme der zweiten Zeile sind aber in denen der ersten sämtlich 
enthalten, da die zweite Gleichung in die erste übergeht, wenn man / durch 
l!+1 ersetzt. Demnach fallen beim Wasserstoff alle Niveaus zusammen, die 
in den Quantenzahlen » und j übereinstimmen: also z.B. 28 mit 2P) 3p. 

2 


| 

fiir die Niv ty=l = — —— (—)E 
” ” ” » J= 2° &=- 2 (=) E. 


mit 3d, usw. 


T 
Als korrigierte Gesamtenergie des Wasserstoffatoms findet man, wenn 
das äußere Feld gleich Null ist: 


Das ist aber nichts anderes als eine Entwicklung der exakten Sommer- BE = 
feldschen Feinstrukturformel nach « Z/n bis zu Gliedern von vierter Ordnung. N oe 
Im Magnetfelde tritt die Alkalinatur des Wasserstoffspektrums klar hervor. Mare 


Die zufällige Form der Koeffizienten A und @ hat jetzt für die Gl. (3,30) keine Re 


Folgen mehr. Die Komponenten der zusammenfallenden Terme, wie z. B. 25, mat ae 

2Pp, verhalten sich durchaus selbständig: Der s-Term spaltet in zwei linear PL \ 

vom Feld abhiingige Niveaus auf, der Py -Term in zwei Terme mit nicht- ap ts, 

linearer Abhängigkeit. Ein Unterschied gegenüber den Alkalien liegt ur n  — a 

dem U all Ni 
gemeinsamen Ursprung aller vier Niveaus. eee 


§ 5. Auswahlregeln und Intensitätsformeln BE? 


Die Eigenfunktionen nullter Näherung hatten wir in Gl. (1,7) in der Form 


(5,1) > a,w" = > + Yy) + 719) 
u 


angesetzt. Die Werte der Koeffizienten x,, y,, sind (bis auf einen gemein- 
samen Faktor) für beliebige Feldstärken aus (3,27) zu entnehmen. Die so be- 
stimmten Funktionen erfüllen jedoch zunächst nur die iterierte Diracgleichung. 
Wir müssen deshalb dafür sorgen, daß unsere Lösung auch die lineare Dirac- 
gleichung, welche als die eigentliche Wellengleichung des Elektrons anzu- 
sehen ist, angenähert erfüllt. Um die Entstehung der iterierten Gl. (1,2) aus 
der linearen deutlich hervortreten zu lassen, schreiben wir sie in der Form: 


„Atombau und nme I, 8. 514. 


, 
rad 
h 
vig 
| 
a 
sd 
23 
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ie E-V &, 


> ie E-V k 
| grad — 5 = 0- 


Hierbei stellen die geschweiften Klammern die beiden linearen Dirac- 
operatoren dar. Diese Gleichung besagt einesteils, daß y eine Lösung der 
iterierten Gleichung ist, andernteils aber auch, gr 

(5,3) y= grad — x) wi -= y 
der linearen Gleichung 
E V E 
5,4) grad - + w= 0 
8 he he 
genügt. Zu einer beliebigen Lösung y der iterierten Gleichung erhält man 
also eine Lösung y’ der linearen nach der Relation (5,3), wobei nur noch 
Differentiationen auszuführen sind. 
Da wir in unserem Falle nur eine Näherungslösung zur Verfügung haben, 


dürfen wir uns damit begnügen, den Operator (5,3) durch eine erste Näherung. Er 


zu ersetzen. Dazu beachten wir, daß die E und E, enthaltenden Glieder 
wegen der Größe der Ruhenergie alles andere überwiegen. Wir schreiben also: 


ie. E-V 
| Ir grad — > x) 


N 


und sehen, daß der Operator bis auf den belanglosen Faktor — ze angeniihert | 
durch den Nullteiler 1 + 7, ersetzt werden kann. Die Eigenfunktionen nullter 


Näherung der linearen Diracgleichung sind danach: 


(5,5) 


bn . 


(5,6) 


R = 
| Ss Tu + U Yas Yn) I ’ 


wobei wir 
(5,7) P= (1 +i + 7) 


setzten. Für die Berechnung der Matrixelemente benötigen wir auch die zu 
(5,6) adjungierte Lösung 


(5,8) Waaj =I > + Yu) y". 


In -dieser Schreibweise ist davon Gebrauch gemacht, daß die Koeffi- 
zienten x, y, sich aus dem Gleichungssystem (3,27) als reelle Größen ergeben. 

Mit yw und '/ wollen wir die Eigenfunktionen nullter Näherung be- 
zeichnen, die zwei beliebigen Atomzuständen zuzuordnen sind. In Uberein- 
stimmung mit (5,6) setzen wir dabei: 


: 5 
660 
AA 
(5,2) 
2 
- 
Br. 
>’ 
ed 
h 
| mod 
Ä 
4% 
a 
— 
2 
> 
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[23 


| y= S eu tine yd 


y (X,, + MT. 
73 


yw“ und ¥* sind nur dann verschieden, wenn sich die beiden Zustände a 
in den Quantenzahlen n oder ! unterscheiden. Nach Dirac ist die dem Uber- an 4 


gang y — entsprechende Übergangsdichte: 

Waaj Y= r> (Wy > (X, + Y,) 


Wie man auf Grund der Definition (5,7) leicht einsieht, ist 


Ti Yas r= 0 
und daher vereinfacht sich (5,10) zu 
5,11) Waa 1 v= > X + yp WAT. 
’ adj /4 


Dieser Form der Übergangsdichte entspricht es, daß wir der Koordinate 4 
das Matrixelement 


(5,12) (q) - fa 2 X, + Yu Y,) dt 


zuordnen. 

Nach dieser Darstellung setzt die Berechnung der Matrixelemente und 
damit der Linienintensitäten die Kenntnis der Funktionen wy“, d.h. also die 
Lösung der Schrödingergleichung des Atoms voraus. Die Koeffizienten x,,--- 
werden durch unsere Störungsrechnung geliefert, und zwar sind, wie das 
Gleichungssystem (3,27) zeigt, im Magnetfelde jeweils höchsten zwei von ihnen 
von Null verschieden, nämlich «,, +1 und Ya-l bzw. X,,_, und Y,,, -1)- 


Daher ergibt die Ausführung der in (5,12) vorkommenden Summationen: 


1 , 1 


(5,13) — 
+ Ym Yu ye dt. 


Für g ist der Reihe nach 


rsinde®?, 
V2 


z=rcos# 
einzusetzen. Die Integrale (5,13) sind, was ihre Winkelanteile betrifft, die- 
selben, wie sie bei der Berechnung des Rotationsspektrums zweiatomiger Mole- 
küle auftreten‘). Es gelten daher wie dort die Auswahlregeln für die Quanten- 
zahlen m und I: 


V=l+1, m=m+1 oder m=m. 


1) Vgl. A. Sommerfeld, ,,Wellenmech. Ergänzungsbd.“ S. 62ff. f 


Er 
4 
~ 
| 
i 
| 
| 


Während die Übergänge m’= m +1 zu den rechts bzw. links zirkular 
polarisierten o-Komponenten des Zeemaneffekts Veranlassung geben, sind die 
Übergänge m’= m für die linear polarisierten 1-Komponenten verantwortlich. 

Wir bemerken ausdrücklich, daB eine entsprechende Auswahlregel für 
die innere Quantenzahl j in unseren Rechnungen nicht enthalten ist. Tat- 
sächlich ist die Bedingung 
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(5,15) 


nur in sehr starken und sehr schwachen Feldern erfüllt. Die Formeln (5,13) 
geben daher auch Rechenschaft von den bei mittleren Feldstärken auftretenden 
„verbotenen Linien“, wie sie von Paschen und Back’) in der I. Nebenserie 
des Funkenspektrums von Calcium entsprechend einem Übergang p, > d, 
beobachtet wurden. 

Zur expliziten Berechnung der Intensitäten aus den Matrixelementen (5,13) 


muß man die Eigenfunktionen wy“ auch in ihren Radialteilen kennen. Diese 
ergeben den gemeinsamen Faktor 
8 8 
oo 
nl nl 
0 


Da F nur von n, 1 und n’, I’ abhängt, hat es denselben Wert für zwei Linien, 
bei denen sich die Ausgangsniveaus ebenso wie die Endniveaus nur in j und m 
unterscheiden, d. h. für solche Linien, die im feldlosen Zustand nur durch 
Feinstruktur getrennt sind (wie z. B. alle 10 Zeemankomponenten der beiden 
gelben Na-Linien). Für die relative Intensität solcher Linien ist also F be- 
langlos. 

Da die Gl. (3,27) nur die Verhältnisse «,, 4/Ym bzw. X, +4/ 


festlegen, müssen wir zur Bestimmung der absoluten Größe der Koeffizienten 
noch die Normierung der Aufenthaltswahrscheinlichkeit heranzuziehen. Setzen 
wir in (5,11) #=w, so erhalten wir die zum Zustand w gehörige Wahr- 
scheinlichkeitsdichte: 


Da die gesamte Aufenthaltswahrscheinlichkeit für ein im Zustand y be- 
findliches Elektron gleich 1 ist, so ergibt sich daraus durch räumliche Inte- 
gration mit Rücksicht auf die Normierung (1,6) der Grundlösungen: 


= 


Die Berechnung der Intensitäten geht nunmehr folgendermaßen vor sich: 
Zu einer gegebenen Feldstärke H sind aus (3,30) die Energiestörungen & und €’ 
des Ausgangs- bzw. Endniveaus zu entnehmen [was allerdings die Kenntnis der 
Dublettabstände (21 + 1) Q,, bzw. (2/’+ 1) Q,,, voraussetzt]. Aus (3,27) er- 
hält man dazu die Koeffizienten Im+t Ym — 4 und X,, +4 Yin’ 1) die 


(5,17) 


dann nach (5,17) zu normieren sind. Damit können die Matrixelemente (5,13) 


1) 


tt 
g 
‘wy 
Al 
| 
k we 
| 
| 
22 
2 
= 1 +4 y™ rs + Yn 2y 
"3 2 
2 
= 
: 
F. Paschen u. E. Back, Physica 1921, S. 
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berechnet werden. Das Quadrat des absoluten Betrags eines Matrixelements 
ist der bei dem betreffenden Übergang ausgestrahlten Energie proportional. 
Die dabei zu beobachteten Linienintensitäten sind dadurch in ihren Verhält- 
nissen ebenfalls bestimmt, aber nach den Gesetzen der Dipolstrahlung je nach 
der Beobachtungsrichtung verschieden '). 

Eine Berechnung der Intensitäten im Zeemaneffekt für beliebige magneti- 
sche Feldstärken ist bereits in der auf S. 658 erwähnten Arbeit von Sommer- 
feld enthalten. Zu den gleichen, entsprechend verallgemeinerten Resultaten 
gelangten dann unter Zuhilfenahme eines Spinmodells Heisenberg und 
Jordan nach quantenmechanischen, Darwin nach wellenmechanischen Me- 
thoden ?). 

Bei schwachen Feldern liefern die Gl. (5,13) die sehr einfachen Inten- 
sitätsformeln, die Hén1*) in diesem Grenzfall für beliebige Multiplizitäten ab- 
geleitet hat. Unter „schwachen Feldern“ sind dabei solche zu verstehen, bei 
denen zwar schon eine Zeemanaufspaltung, aber noch keine Beeinflussung der 
Übergangswahrscheinlichkeiten durch das Magnetfeld wahrnehmbar ist. Wir 
müssen also die Koeffizienten x, --- für H = 0 berechnen, um zu den Hönl- 
schen Formeln zu gelangen. Die Koeffizienten werden dabei von der Größe Q 
| unabhängig. Die gleiche Vereinfachung ergibt sich in dem anderen Grenzfall 

des Paschen-Backeffekts. Er ist dadurch gekennzeichnet, daß in (3,27) Q 
gegen @ zu vernachlässigen ist. Auch hier gehen also die individuellen Kon- 
stanten Q der einzelnen Dubletts in die Intensitätsverhältnisse nicht ein. 


Zusammenfassung 

Zur Berechnung der Feinstruktur und des Zeemaneffekts der 
Alkalien wird nach einem von Sommerfeld und Maue angegebenen 
Verfahren die iterierte Diracgleichung näherungsweise integriert. 
Die Eigenwertstörungen werden in zweiter Näherung nach der 
Schrödingerschen Störungstheorie ermittelt. Es ergeben sich 
dabei die Voigt-Sommerfeldschen Formeln die Zeemanaufspal- 
tungen von Dublettsystemen. Mit Hilfe der Eigenfunktionen nullter 
Näherung werden die Auswahl- und Intensitätsregeln im Zeeman- 
effekt diskutiert. 


Meinem hochverehrten Lehrer, Herrn Prof. A. Sommerfeld, 
möchte ich für die Anregung zu dieser Arbeit und für seine immer- 
währende Unterstützung bei ihrer Durchführung meinen tiefsten 
Dank aussprechen, ebenso auch Herrn Dr. habil. A. W. Maue für 
viele aufschlußreiche Diskussionen. 


2) Vgl. die Zitate auf S. 643. 
3) F. Hönl, Ztschr. f. Phys. 31. 8. 340. 1925. 


Höllriegelskreuth b. München, Ges. f. Lindes Eismaschinen- 
A.-G. x 


(Eingegangen 30. September 1936) 
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Logarithmische und lineare Lautstärkenskala Pr 
Von W.Bürck, P.Kotowski und H.Lichte 


Vor zehn Jahren hat Barkhausen!) die Lautstärkeeinheiten des 
Wiens und des Phons vorgeschlagen, wobei die erste Einheit einer 
linearen und die zweite einer logarithmischen Skala zugeordnet 
gewesen ist. Die Begründung der logarithmischen Skala war dabei 
der große Umfang der Lautstärkemaßzahlen bei linearer Skalen- 
teilung, der durch das logarithmische Maßsystem erheblich ver- 
kleinert werden konnte. Diese zweite Einheit ist wiederholt mit 
der Annahme des Weber-Fechnerschen Gesetzes bzw. seiner 
strengen Gültigkeit in gewissen Bereichen verbunden worden?), so 
daß heute mit den Lautstärkemaßzahlen der logarithmischen Skala 
von vielen Leuten unklare und verschwommene Vorstellungen über 
die „empfundene Lautstärke“ verbunden werden. Diese Annahme, 
daß die logarithmische Skala dem Ohr besonders entspräche, hat 
anscheinend mit dazu geführt, daß sie in letzter Zeit fast aus- 
schließlich benutzt wird, während die von Barkhausen vorge- 
schlagene, dem normalen Gebrauch näherliegende lineare Skala fast 
vergessen ist. Auch Wien?) bediente sich bei seinen grundlegenden 
Lautstärkemessungen einer linearen Skala. Er wählt den logarith- 
mischen Maßstab nur, um eine graphische Darstellung zu ermög- 
lichen, weist allerdings auch auf die Berechtigung dieser Einteilung 
für den Fall hin, daß das Weber-Fechnersche Gesetz Gültigkeit 
besitzt. Die Lautstärkeempfindung ist aber wie jede Empfindung mit 
physikalischen Methoden grundsätzlich nicht meßbar, so daß man die 
„richtige Skalenteilung‘‘ mit Hilfe der Empfindung nicht prüfen kann. 
Mit dem Ohr messen kann man lediglich die Gleichheit zweier 
Lautstärken, genau wie man z. B. in der Optik mit dem Auge nur 
die Gleichheit zweier Helligkeiten feststellen kann, so daß man eine 
willkürliche Lautstärkeskala für Schallereignisse einer bestimmten 
Art z. B. den reinen Sinuston von 1000 Hz festsetzen muß, worauf 
man zu einem unbekannten Schallereignis den gleichlauten Ton der 


1) H. Barkhausen, Ztschr. f. techn. Phys. 7. 8.59. 926. = =| 
2) Vgl. z.B. R. Berger, Schalltechnik. 6. S. 17. 1933. 
3) M. Wien, Phys. Ztschr. 4. S. 69. 1902. 
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Normalskala bestimmt und dessen willkürlich definierte „Lautstärke“ 
als Lautstärke des unbekannten Schallereignisses übernimmt. Da 
man die Skala willkürlich teilen kann, wird man zweckmäßig ein- 
fache Gesetze zur Einteilung wählen, wobei die lineare Skala die 
einfachste ist, die außerdem jedem, auch dem physikalisch und 
mathematisch Ungebildeten leicht verständlich ist. Hält man den 
Vorzug eines engen Zahlenbereichs der möglichen Lautstärken für 
wichtig genug, so kann man die logarithmische Skala wählen. Man 
führt damit aber etwas sogar dem Physiker Ungewohntes ein und 
verzichtet auf ein Verständnis der akustischen Maße durch Nicht- 
fachleute. Insbesondere sind die gewöhnlichen Rechenoperationen 
der Addition und Multiplikation nicht mehr in der gebräuchlichen 
Art ausführbar. Es wird beispielsweise einem Nichtfachmann 
schwer erklärlich sein, daß die Aufstellung von zwei gleichlauten 
Schallquellen einer Steigerung der Lautstärke gegenüber einer Schall- 
quelle um einen konstanten Phonbetrag unabhängig von der Ab- 
solutlautstärke gleichkommt. ® 

Die Verhältnisse in der Akustik sind denen in der Optik weit- = 
gehend gleich. Obwohl das Auge nach Ansicht vieler Leute log- 
arithmisch registriert, hat kein Physiker und Beleuchtungstechniker == 
eine logarithmische Skala für Helligkeiten benutzt. Es hat sich hier ee 
vielmehr trotz des ähnlich großen Meßzahlenumfangs wie bei Schall- 
messungen nur die einfache lineare Skala durchgesetzt und bestens 
bewährt. Daneben wird noch eine logarithmische Skala für einzelne 
Fälle (Sterngrößen in der Astronomie, Schwärzung in der Photo- 
graphie) benutzt. 

Eine Skala nach der Empfindung aufzubauen, wird niemals ein- 
deutig gelingen und ist außerdem für den praktischen Gebrauch 
unzweckmäßig. Will man aber die „empfundene Lautstärke“ oder 
die „Lautheit“ als Einteilungsgrundlage in der Akustik wählen, so 
ist man gezwungen, willkürliche Annahmen über den Gang der 
Lautheit mit der Schallintensität zu machen. Fletcher’), der diese 
Frage sehr eingehend bearbeitete und auch die Ergebnisse; fremder 
Untersuchungen zusammenfaßte, kommt zu folgendem Ergebnis: 
Drei verschiedene Methoden führen zu annähernd derselben Laut- 
heitskala. Erste Untersuchungsmethode mit folgender Annahme: 
Bei gleicher Empfindlichkeit beider Ohren ist die bei einohrigem 
Hören empfundene Lautstärke gleich der Hälfte der bei zweiohrigem 
Hören empfundenen. Zweite Untersuchungsmethode mit der An- 
nahme: Zwei Töne, die jeder für sich die gleiche e 


Re 


1) H. Fletcher, Journ. Frankl. Inst. 22. 
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Lautstärke hervorrufen, ergeben zusammen die doppelte Lautstärke, 


falls sie so gewählt werden, daß keine Verdeckungseffekte auftreten. 2 B 


Dritte Untersuchungsmethode: Ohne physikalische Annahme wird 
bei verschiedenen Ausgangslautstärken die subjektiv empfundene 
Lautstärkeverdoppelung registriert. Die beiden ersten Methoden 
lassen sich auf die gemeinsame Grundannahme zurückführen, dab 
die Anzahl der dem Gehirn zugeleiteten Nervenimpulse proportional 
dem zugeführten Schalldruck ist, und daß, zum mindesten im Be- 
reich 1:2, die Empfindung linear von der Impulsanzahl abhängt. 
Die Abhängigkeit der Empfindung von der Impulszahl kann aber 
mit physikalischen Mitteln nicht nachgeprüft werden. Die dritte 
Untersuchungsmethode ist nur möglich, wenn sich die Meßpersonen, 
durch längere Übung befähigt, einen Empfindungsmaßstab zurecht- 
gelegt haben, der sich jedoch jeder physikalischen Kontrolle entzieht. 

Daß der Empfindungsmaßstab tatsächlich ein Geschmacksurteil 
darstellt, ist aus der Tatsache ersichtlich, daß z. B. ein Schall- 
ereignis A leiser als ein anderer Schallvorgang B erscheint, dieser 
Vorgang B wieder leiser als C, trotzdem aber in reproduzierbarer 
Weise von manchen Beobachtern A lauter als C empfunden werden 
kann}, Steudel folgert daraus, daß der Lautstärkemaßstab dieser — 
Gruppe von Beobachtern falsch sei; nach seiner Meinung müßten 
diese Meßpersonen sich entweder einen anderen Lautstärkemaßstab — 
angewöhnen, oder sie kämen für akustische Messungen als Beobachter 
nicht in Frage. Die Ansicht, daß ein Empfindungsmaßstab falsch 
sei, ist aber ebenso unberechtigt wie die, man könne einen anderen 
Maßstab als „richtig“ bezeichnen, da nicht einzusehen ist, warum 
für Empfindungen die Regeln der Algebra Gültigkeit besitzen sollen. 

Man muß also irgendwelche willkürlichen Annahmen bei der 
Aufstellung einer Empfindungsskala machen. 

Unter den Voraussetzungen, die Fletcher gemacht hat, ge- | 
langt man zu einer ,,Lautheitsfunktion“, die je nach der Art der 
Schallvorgänge etwas verschieden ausfällt, im allgemeinen aber den 
in der Abbildung als Kurve a dargestellten Verlauf hat. In der 
Abbildung ist gleichzeitig eine Linearteilung (Kurve b) und eine 
logarithmische Abhängigkeit (Kurve c) so eingetragen, daß innerhalb 
der . dargestellten 5 Zehnerpotenzen des Reizes (Schalldruck) der 
größte Fehler in der Ordinatenrichtung ein Minimum wird. Dieser 
größte auftretende Fehler ist bei linearer Skala etwa der Faktor 3,5, 
bei logarithmischer Skala etwa der Faktor 20. Die Fletchersche 
Lautheitskurve schmiegt sich also wesentlich besser der linearen 


1) U. Steudel, u. 41. 8. 116. 
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als der logarithmischen Teilung an. Je größer der Bereich gewählt 
wird, um so fehlerhafter wird die logarithmische Teilung; in ge- 
nügend kleinen Bereichen läßt sich ja jede Teilung mit der Laut- 
heitskurve in Einklang bringen. Selbst unter der Voraussetzung 
also, daß man die physikalische Skala der Lautstärke nach der 


Lautheit einteilen wollte, hat die Linearteilung vor der logarith- 
mischen den Vorrang, da das Weber-Fechnersche Gesetz die 
empfundene Lautheit schlechter wiedergibt als die einfache An- 
nahme eines linearen Zusammenhangs. 

Die oft für die Gültigkeit des Weber-Fechnerschen Gesetzes 
zitierte Erscheinung, daß z. B. bei einer Reihe von gleich großen, 
nach Helligkeit gestuften Flächen der Helligkeitszuwachs proportional 
zur vorher vorhandenen Helligkeit gewählt wird, sagt nur aus, daß 
die Steigerung eines Sinneseindruckes nach einer geometrischen 
Reihe eine ästhetisch befriedigendere Wirkung auf den Durchschnitts- 
menschen ausübt, ihm als „gleichmäßiger“ erscheint, als eine Stei- 
gerung nach einer arithmetischen Reihe. Hieraus einen Zusammen- 
hang mit der Empfindungsskala herzuleiten, ist willkürlich. 

Die logarithmische Lautstärkeskala hat auch deswegen die 
lineare Skala verdrängt, weil zur Zeit Schwachstromtechnik und 
Akustik fast überall gemeinsam behandelt werden und die Schwach- 
stromtechniker an logarithmische Maßstäbe von der Kabel- und Ver- 
stärkertechnik her gewöhnt sind, in der derartige Maßstäbe viel 
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Fällen zweckmäßiger sind als lineare. Dabei wird häufig übersehen, 
daß die in der Schwachstromtechnik üblichen logarithmischen Maße 
neben häufig gebrauchten linearen bestehen, so daß man sich dort 


Skala bereits ungebräuchlich geworden. Es erscheint also mindestens 
der Gebrauch einer linearen Skala neben der logarithmischen als 
berechtigt, und dies um so mehr, als von den bisher vorgeschlagenen 
logarithmischen Skalen bereits zwei, nämlich die Skala der alten 
Phons und der S.U. (sensation units) nach kurzem Gebrauch nicht 
befriedigten. Erwägenswert scheint in diesem Zusammenhang der 
Vorschlag von Davis!), für alle Zweige der Physik eine neue ein- 
heitliche logarithmische Skala einzuführen, die durch besondere Bei- 
wörter auf die verschiedenen Arbeitsgebiete bezogen werden könnte. 


1) A.H. Davis, Proc. Phys. Soc. London. 46. S. 631. 1934. 
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